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Abstract - This paper continues the study, initiated in |28l of the validity of the 
Zakharov model describing Langmuir turbulence. We give an existence theorem for a class 
of singular quasilinear equations. This theorem is valid for well-prepared initial data. We 
apply this result to the Euler-Maxwell equations describing laser-plasma interactions, to 
' obtain, in a high-frequency limit, an asymptotic estimate that describes solutions of the 

Euler-Maxwell equations in terms of WKB approximate solutions which leading terms are 
, solutions of the Zakharov equations. Because of transparency properties of the Euler- 

f"*^ ' Maxwell equations put in evidence in |28|. this study is led in a supercritical (highly 

nonlinear) regime. In such a regime, resonances between plasma waves, electromagnetric 
waves and acoustic waves could create instabilities in small time. The key of this work is 
i^h | the control of these resonances. The proof involves the techniques of geometric optics of 

Joly, Metivier and Rauch |13II14| . recent results of Lannes on norms of pseudodifferential 
^ ■ operators |I]>., and a semiclassical, paradifferential calculus. 
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1 Introduction 



We describe solutions of initial value problems for quasilinear, hyperbolic 
systems of the form, 

dtu + -^rA(e,eu,ed x )u = -B(u,u), (1.1) 

e z e 

in the high-frequency limit e — > 0. 

In (jl.lj) . A is a symmetric hyperbolic, differential, or pseudo-differential 
operator; the singular source term B is bilinear. The unknown u £ has values 
in W 1 . It depends on time t G R + and space x G and is subject to the 
initial condition, 

u £ (0,x) = a £ (x), (1.2) 

where a £ is a bounded family in H s (JH d ), for some large s. 

In this setting, the existence, uniqueness and regularity of solutions to 
(jTTTjl - lfOj) for fixed e > is classical. 

The limit e — > is singular in two ways: first, solutions develop fast 
oscillations in time, with frequencies of typical size 0(l/e 2 ); second, the 
amplitude O(l) of the initial datum is large, hence the singular source term 
B/e could create instabilities in small time 0{e). 

Under appropriate assumptions, we prove the existence of solutions to 
(|1.1|) - ([1.2|) over time intervals independent of e, and their stability under 
initial perturbations of the form e k °ip £ , where ko is large enough, and <*p £ is 
bounded in a semiclassical Sobolev space, in particular, may contain fast 
oscillations of the form e tkx f £ . 

We show that our assumptions are satisfied by the Euler-Maxwell equa- 
tions describing laser-plasma interactions. This implies in particular that, 
in a high-frequency limit, solutions of the Euler-Maxwell equations are well 
approximated by WKB approximate solutions which leading terms are so- 
lutions of the Zakharov equations. 

Our assumptions and results are precisely stated in section 12.41 

1.1 Discussion: weakly nonlinear versus highly nonlinear ge- 
ometric optics 

Consider 1)1.1(1 . and assume for instance that A is a differential operator of 
the form 

A(e,£u,££) = A Q {e) +eAi(e,g) +e 2 A 2 (e,u^), (1.3) 
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where Ai and A2 are linear in £, and Ao, A\ and A2 are hermitian, for 
all £,u,£. Suppose that the family of initial data has the form u £ (0,x) = 
e p a £ (x), where a 6 is a bounded family in H s (M. d ). 

The weakly nonlinear regime corresponds to p = 1. In this regime, the 
classical existence proof provides a maximal existence time t*(e) that satis- 
fies liminf e ^o t*(e) > 0. Indeed, if one lets v £ = eu £ , then the initial datum 
for v £ is O(l), and the equation in v £ is 

d t v £ + i?A(e,ev £ ,£d x )v £ = B(v £ ,v £ ). 

The classical H s energy estimate for quasilinear symmetric hyperbolic op- 
erators then yields 

\v £ {t)\ H s < \a £ \ H s +C f \v £ (t')\ H sdt', 
Jo 

where C depends on |u £ |^n,oo, then, with Gronwall's lemma, the uniform 
bound 

\v £ {t)\ H ° < \a £ \ H se ct . 

This uniform estimate is the key of the proof of the existence of a solution 
over a time interval independent of e. 

On the contrary, whenp = 0, the initial data have a large amplitude 0(1). 
In this regime, the maximal existence time a priori satisfies t*(e) = 0(e), 
and in particular lim inf £ _>o t* ( £ ) = 0. The H s energy estimate gives indeed 

|u e (t)|ff' < \a £ \ H * + - I \u £ (t')\ Hs dt', 
£ Jo 

whence, 

\u £ {t)\ H s < \a £ \ H se ct/£ . 

This shows that the singular term B/e may cause the solution to blow- 
up in small time. This highly nonlinear (or supercritical) regime is called 
oscillations fortes, or strong oscillations, by Cheverry, Gues and Metivier in 

Our goal is to state conditions on (|1.1|) - (|1.2|) that are sufficient to have 
liminf e _>o t*(e) > in the regime p = 0, and that are satisfied by the Euler- 
Maxwell equations. 
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1.2 The Euler-Maxwell equations and the Zakharov equa- 
tions 

This paper is in the direct continuation of |281 I8j. The underlying physical 
context is the study of laser-plasma interactions; in particular, the question 
of the rigorous derivation of the Zakharov model from fundamental equa- 
tions. 

We take here as a system of fundamental equations the Maxwell equa- 
tions coupled with the Euler equations [23] > 

d t B b + cV x E b = 0, 

8 t E b - cVx# = 4vre((n + n\)v\ - (n + n b )v b )), 

m e (n + n b e )(d t v b e + (v b e ■ V)v b e ) = - le T e Vn b e - e(n + n b e ){E b + -v\ x B b ), 

c 

mi{n + n b )(d t v b + (y\ ■ V)v\) = - 7i TiVnJ + e(n + n b e ){E b + -v\ x B b ), 

c 

d t n\ + V • ((n + n\)v\) = 0, 
d t n\ + V • ((n + n b e )v b ) = 0. 

The variables are B b , E b the electromagnetic field, v b e ,v b the velocities of the 
electrons and of the ions, and n e and n b the density fluctuations from the 
equilibrium no- The first two equations are Maxwell's equations describing 
the time evolution of the electromagnetic field, the next two equations are 
the equations of conservation of momentum for the electrons and the ions, 
and the last two equations are the equations of conservation of mass for the 
electrons and the ions. The electric charge of the electrons is — e; to simplify, 
we assume that the charge of the ions is +e. The parameters are m e and rrii 
the masses of both species, j e and 7, the specific heat ratios of both species, 
Tj and T e the temperatures of both species and no the (assumed constant 
and isotropic) density of the plasma at equilibrium. 

In the above system, Maxwell's equations are coupled to Euler's equa- 
tions by the current density term in the right hand side of the Ampere 
equation and by the Lorentz force in the right hand side of the equations of 
conservation of momentum. 

The additional divergence equations 

V-5 b = 0, V ■ E b =4Tre(n b e -n b ), (1.4) 

are satisfied at all times if they are satisfied by the initial data. A brief 
discussion of the relevance of this model is given in [^Sj. We work in this 
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paper on the non-dimensional form of these equations introduced in j^H] : 

d t B + X7 x E = 0, 

d t E - V x B = ~(1 + nl)v e - ~(1 + 4) Vi , 

£ £ Up 



(EM)' < 



9fu e + 8 e (v e ■ V)v e = -9 e 



Vn 



1 



~(E + 6 e v e xB) 



l + n e £ 
dtnl + e V • ((1 + 4)v e ) = 0, 



1 + n\ £Ve 

dt4 + OiV ■ ((1 + njj)^) = 0. 
In (EM) tt , the variable is 

J = (B,E,v e ,4,Vi,4) el 14 . 
The change of variables for small amplitudes, 



1 + 4 = e n °, l + nf = e"\ 



(1.5) 



leads to the system, 



(EM) < 



d t B + V x E = 0, 
d t E-VxB = -(e n "v e -^e ni v i ), 

£ V e 

d t v e + 8 e (v e ■ V)v e = -6 e Vn e - ~(E + 6 e v e x B), 
d t n e + 6> e V • v e + 6 e (v e ■ V)n e = 0, 

dm + Qiivi ■ V)vi = -a 2 uiVni + -±-(E + 6m x B) 



£0, 



d t rii + OiV ■ Vi + OiVi ■ Vrii = 0. 
In (EM), the variable is 



n ; 



u = (B,E,v e ,n e ,Vi, — ) G 
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where (B,E) £ M 3+3 is the electromagnetic field, (v e ,Vi) £ M 3+3 are the 
velocities of the electrons and of the ions and (n e , rtj) £ are the fluctu- 
ations of densities of both species. The variable u depends on time t € M+ 
and space x£l 3 . The small parameter e is 

1 



£ :- 



LO pe to ' 
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where uj pe is the electronic plasma frequency: 



47Fe2no (1.6) 



m e 



and to is the duration of the laser pulse. A typical value for e in realistic 
physical applications is e ~ 10~ 6 . The parameters a,9 e and 9i are 



c\ m e c V mj T e 

Typically, # e ~ 10~ 3 . Because the ions are much heavier than the electrons, 
9{ and a are much smaller than 9 e . We consider the specific regime 

9i = e, (1.7) 

and we look for solutions to (EM) with initial data of size 0(e) defined over 
diffractive times 0(1/ e), that is, we make the ansatz 

u(t, x) := eu(et, x). (1.8) 

Written as a system of equations in the variable it, (EM) takes the form 

The Zakharov system is a simplified model for the description of the 
nonlinear interactions between E, the envelope of the electric field, and n, 
the mean mode of the ionic fluctuations of density in the plasma, 

-- nE, 

This model was derived by Zakharov and his collaborators in the seven- 
ties [IHj- It describes nonlinear interactions between high-frequency, electro- 
magnetic waves and low- frequency, acoustic waves. In (Z), the Schrodinger 
operator is the classical three-scale approximation of Maxwell's equations 
jl.Sj ; the wave operator is the classical long- wave approximation of the Euler 
equations. The nonlinear term in the right-hand side of the equation in E 
directly comes from the current density term in the Ampere equation. The 
term A|.E| 2 comes from the convective terms and the nonlinear force term 
in the equations of conservation of momentum. A WKB expansion of the 
(EM)" system is performed in section r3.ll In the nonlinear regime of our 
interest, the limit system is (Z). 
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1.3 Description of the results 

We extend here the results of [281 15]. as we show that, in the high-frequency 
limit and in a highly nonlinear regime, solutions of the Euler-Maxwell equa- 
tions are well-approximated by solutions of the Zakharov equations. In 
section |21 under appropriate assumptions, we prove the following theorem: 

Theorem 1.1. The unique solution to (|1.1|) - (|1.2[) is defined over a time 
interval independent of e. 

Then we apply this result to the Euler-Maxwell equations, to obtain: 

Theorem 1.2. In the high-frequency limit, solutions of the Euler-Maxwell 
equations initiating from polarized, large- amplitude initial data, are well ap- 
proximated by solutions of the Zakharov equations initiating from nearby 
initial data, in the sense that there exists to > and C > 0, independent of 
e, such that, 

sup sup(|£ - (Ee iuJ " et/£2 + c.c.)| + \n £ - en\) < Ce 2 , 

0<t<t x 

for e small enough, where E e and n £ represent the electric field and the ( elec- 
tronic or ionic) fluctuation of density in the solution of the (EM) system, 
and E,h, is the solution of (Z). 

We show in sectionOJhow Theorem ll.2l follows from Theorem ll.il Precise 
statements are given in sect ions 12 . 41 ( Theror em . 1 5 1) and 13 . 21 (Theorem 13 . 1 1) . 

1.4 Outline of the proof 

The proof of Theorem 11.11 ( section 12. 5 J) goes along the following lines: the 
construction of a precise, regular, polarized approximate solution defined 
over a time interval independent of e, the preparation of the system, and 
finally, the control of the resonant interactions of oscillating waves. 

1.4.1 Existence of an approximate solution 

In section El we consider the initial value problem for the (EM) system, and 
show that it takes the form 1(1 - 1)1.21) . We construct a WKB approximate 
solution, 

u e a = u + eui + . . . , (1.9) 
under an assumption of polarization for the initial data. 
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The polarization condition is a well-preparedness condition; it is neces- 
sary for the construction of a WKB approximate solution that is continuous 
in time, at t = 0. Consider indeed u e a in the form 1)1. 9|) . a solution to (Jl.lj) - 
(|1.2|) . where the operator ^4 has the form l)1.3|) . with *4o(0) 7^ 0, as in the 
(EM) system. The limits e — > yields 

A (o)«o = 0, 

then, if ito is continuous at i = 0, 

A(0)a° = 0. 

The above condition is the well-preparedness assumption for the initial da- 
tum a £ . For the (EM) system, it takes the form 

a = (0 K 3 , E + (E)*, —E — —E*, Or, r3 , Or) , 

for some fundamental frequency u>, defined in terms of Lo pe , and some com- 
plex amplitude E (above, (E)* denotes the complex conjugate of E). 

WKB solutions to the Euler-Maxwell equations, initating from highly- 
oscillating, and well-prepared, initial data, are considered in [2Sj. It is shown 
in [1H] that, 

1) the (EM) equations satisfy transparency properties, that is, null condi- 
tions for coefficients describing constructive interactions of character- 
istic waves. As a result, the weakly nonlinear (in the sense of section 
11.1(1 approximation of the (EM) system is a linear transport equation, 
and 

2) WKB solutions of the (EM) equations, initiating from large-amplitude 
solutions, satisfy, in the high-frequency limit e — > 0, systems of the 
form 

J i{8 t + cd z )E + A ± E = nE, 
( M (d?-A ± )n = A ± (\Ef), 

where z is the direction of propagation of the laser pulse, and is 
the Laplace operator in the transverse directions. 

The approximate solution that is constructed in section El satisfies the 
ansatz, 

u e a (t,x) = U e a (t lXj %). (1.10) 
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In particular, there are three times scales. This is consistent with the 
well-known fact that the Schrodinger equation is an approximation of the 
Maxwell equations in the diffractive limit (that is, t = 0(1) and oscillations 
with frequencies 0(l/e 2 )). Note, however, that the wave equation, also 
present in the (Z) system, is an approximation of the Euler equation in 
the geometric optics limit (that is, t = O(l) and oscillations in 0(l/e)). 
The third scale is actually built in the Euler equations by the "cold ions" 
assumption 9% = e. 

In (jl.lOj) . the profiles are purely time-oscillating. In particular, the initial 
data do not have fast oscillations. The limit system is (Z) , that is, a 
Zakharov system with zero group velocity (see the characteristic variety 
pictured on figure |2J). Such waves are called plasma waves in the physical 
literature. The approximate solution has the form 

< = (n ,ie^ /£2 + <ie-^ 2 ) + e(u lfi + ...)+ e 2 (, . . ), (1.11) 

where the components of no 1 and u± o satisfy (Z). 

For general equations of the form (|1.1|) . studied in section |2J the existence 
of an approximate solution is assumed (Assumption I2.l1j) . 

1.4.2 Preparation of the system 

Given a precise approximate solution u £ a of the form (|1,11|) , we look in section 
I2.5.1l for the exact solution u £ as a perturbation of u £ a , 

U 6 = U E a + £ k U S . 

The initial condition is u £ (t = 0) = a £ + e k °ip £ , where (p £ has a high Sobolev 
regularity and ko is large enough. In the definition of u £ , k is chosen in terms 
of ko. One assumes that u £ a is accurate at an order Iq, much larger that k$. 
The equation in ii £ has the form, 

d t u E + ^A{e{u e a + e k ii e ))u £ = ^(B(u £ a ,u e ) + B(u £ ,u £ a )), (1.12) 

The propagator A/e 2 is hyperbolic fAssumption l2.10"|) . thus generates highly 
oscillating waves, with frequencies of typical size 0(l/e 2 ). We write the 
spectral decomposition of A as follows: 

A = J ^2 ^kg n kg + ^acll a c. 
kg ac 

The real eigenvalues A^g are called Klein-Gordon modes, while the real eigen- 
values A ac are called acoustic modes. The characteristic variety for the (EM) 
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system (that is, the union of the graphs £ i— ► Aj(£), at u = 0, for j = kg and 
j = ac) is pictured on figured In particular, A^g ~ 0(1)£, while A ac ~ s£, 
a consequence of the cold ions hypothesis ()1.7|) . 

The Klein-Gordon waves generated by A interact with the highly os- 
cillating approximate solution, through the convection term, and through 
the source term B. These interactions create low-frequency waves, which 
can be seen as source terms in the equations for the components of the so- 
lutions in the directions of the acoustic modes. Thus low-frequency, and 
high-frequency waves are propagated. The Zakharov system pretends to 
describe how these waves interact. 

Equation (|1.12|) . together with an initial datum of size 0(1), can be 
likened to an ordinary differential equation, 

. ia 1 2 

y +-ny = -y , 

e z e 

with an initial datum y(0) = y$. The singular source term in the right-hand 
side may cause the solution to blow-up in small time 0(e), but exponential 
cancellations are expected to happen because of the rapid oscillations created 
by the term in 1/e 2 . 

To investigate these exponential cancellations, it is natural to project the 
source term B/e over the eigendirections of A. Let the total eigenprojectors, 



We compute 



n = ^n kg , n s = ^n ac . 

kg ac 



n sn ii bu s \ = ( * o(i) 
u s bu u s bu s ) \ ** * 



Because constructive interaction of waves between low- and high-frequency 
do occur for the (EM) system (see figure^), the 0(1) term in the right 
block of the above interaction matrix can be interpreted as an absence of 
transparency. 

We then rescale the solution (section 12.5. 2|) as follows, 

v £ ■= (n ii £ ,-u s u £ ). 



e 



The equation in v £ has the form 



d t v £ + ^Av £ = ^Bv £ + -Dv £ + 0(l)v £ 
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with the notations, 

"-en). °-{i°.y 

In the equation for v, the propagator A is diagonal and the leading source 
term B is nilpotent. The system is prepared. 



1.4.3 Control of the constructive interactions of waves 

In a third step (sections 12.5.31 to 12. 5. 5j) . the singularity in the right-hand side 
in v e is reduced. Consider a change of variable in the form 



w £ := (Id + iV)-V, N 




N 



The equation satisfied by w £ is 



d t w £ + ^Aw £ = 4j (-B — [e 2 d t + A, jV])u? s + ^Dw £ . 

We look for N solution of the homological equation, 

B-[e 2 d t + A,N] = 0(e 2 ), 

This equation takes the form = B + 0(e 2 ), where the phase $ is = 
Akg — A ac — oj. The equation $ = is the resonance equation. Its solutions 
are pictured in Figures 01 IU and |SJ The crucial transparency assumption 
(Assumption I2.13|) states that the interaction coefficient B is sufficiently 
small at the resonances, that is, 

\B\<Ce 2 \$\. (1.13) 

The equation in w e becomes 

d t w £ + \aw 6 = -Dw e . 

£ z £ 

A symmetrizability assumption for D ( Assumption 12 eventually allows 
to perform energy estimates (section I2.5.6JI . which yield uniform bounds for 
w £ , and a continuation argument concludes the proof. 

In sectional we describe the resonance equations for the Euler-Maxwell 
equations, and check that an estimate of the form (|1.13|) is satisfied. 
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1.4.4 Technical issues 

The symbols in the spectral decompositions of A are necessarily pseudo- 
differential operators, even when A is differential. They also depend on the 
solution u £ , because the equations are nonlinear. We are naturally led to 
consider pseudo-differential operators of the form 

q{e,x,g) =p(£,v(x),g), 

where v has a Sobolev regularity. The questions of the bounds of the cor- 
responding operators, in Sobolev spaces, and of the existence of a symbolic 
calculus, naturally arise. Lannes recently gave optimal bounds in |15j . For 
an operator of order m, these bounds have the form, 

\\op(p(v,£)u\\H° < C(\v\l°°)(\\v\\h°\\u\\h s o + \\u\\ H s+m). 

A symbolic calculus is available; the operator op{pi(v))op{p2(v)) has the 
symbol 

Pi(v) P 2(v) + Y, d?( Pl (v))d«(p 2 (v)) + ... 

\a\=l 

When the symbols depend on x through the solution u £ , the subprincipal 
symbol depends on d x u £ , and its operator norm in H s depends on ||u e ||#-s+i. 
That is, compositions of such operators lead to losses of derivatives. 

To overcome this difficulty, it is classical to differentiate the equation up 
to order s, and then perform energy estimates in L 2 . 

We now explain why we chose a different approach. 

In the perturbation equations (|1.12jl . all the derivatives are e-derivatives 
(see below). The equation in ed x ii £ (we let d = 1 in this discussion), has a 
singular source term in u £ . The variable U £ := {ii £ , (ed x )u £ ), solves 

<vr + ±msu-)u- = i(£ ° B )v-±( ) u; 

where B is short for B{u £ a ) and B' is short for B{d x u £ a ). Up to 0(e 2 ), the 
commutator in the above equation is ed v A{u £ a ). Following the approach of 
section Tl.4.21 we rescale the solution, by letting 

V £ = (n n e , -U s u £ ,U ed x ii £ , -IL s sd x ii £ ). 

The equation becomes, 

d t V £ + \A(eu £ )V £ = ^BV £ + -DV £ , 
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with the notation, 



B := 



I \ 

u s bu ooo 

o ooo 

V e{n s B' - d v A(u a ))u o n s £>n o / 



Because of the singular rescaling, the commutator H s d v A(u a ))Ho is now a 
singular source term. It cannot be symmetrized in the matrix B_, and might 
well be not transparent. That is, B might not satisfy an estimate of the form 
(|1.13[) . even if B does. In the case of the (EM) system, this term can actually 
be shown to satisfy a transparency estimate of the form (|1.13|) , so that the 
method of differentiating the equation could be applied. To obtain a iJ| 



and notational complications would arise in the normal form reductions of 
sections [2.5. 31 to 12 . 5 . 51 This method would arguably be conceptually simpler 
than the one we chose, namely paradifferential smoothing. 

The paradifferential smoothing of Bony |Sj is another classical way to 
overcome the artificial losses of derivatives that occur in the compositions 
of operators. We denote the paradifferential operators by op 1 ^, where tp is 
an admissible cut-off (see section • It is classical that for an operator of 
order m, 



and the norm of the subprincipal symbol in the composition of two para- 
differential operators depending on is m | u|^/i,oo. 

The setting of our interest is semi-classical, in the sense that the opera- 
tors depend on £ through e£. It is easy to check that the above bounds and 
symbolic calculus can be adapted to this setting (section \2.2\i . In a semi- 
classical setting, subprincipal symbols arising in the compositions come with 
a prefactor e. Because the singularity is in 1/e 2 , this implies that we need 
only consider the principal and the subprincipal symbols. The perturbation 
of the initial data is accordingly assumed to have a semiclassical Sobolev 



estimate, one would have to write a system of size 2n ^ 
with a leading singular source term in a bidiagonal form 




/ B \ 

1 eB' B 

e 2 

\ eB' B J 



op 



(p(v),£)u\\h» < C{\v\ L °°)\ 
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regularity. In particular, it can take the form (p(x)e ikx ' £ , with 99 E H s . The 
final asymptotic estimate (|2.29j) is formulated in H*. It implies an estimate 
in L°°. 

We finally mention a technical point, associated with a lack of regularity 
of the operators involved in the changes of variables described in sections 
11.4.21 and 11.4.31 caused by the fact that the spectral decomposition of the 
hyperbolic operator in the (EM) equations becomes singular for small fre- 
quencies. 

The wave equation in (Z), that comes up as a geometric optics approxi- 
mation of the Euler equations, is associated with symbols in ±e|£|. In par- 
ticular, these symbols are only bounded at the origin. Because resonances 
between Klein-Gordon and acoustic modes occur precisely for small frequen- 
cies |£| ~ e, a smoothing procedure, or the introduction of a cut-off, would 
create large error terms. At this point, we make a crucial use of the fact 
that all the symbols depend on the solution u £ , only through eu £ , and that 
we work on perturbation equations: u e = u £ l + e k ii £ , where k is large enough. 
Because we need to handle symbols only up to 0(e 2 ), we can approximate 
p(su £ ,£^) by p(0, e£) + ed u p(0, s£) • u £ a . This approximate symbol is easier 
to handle, for two reasons. First, it depends on x only through u £ a , the ap- 
proximate solution. Second, it has the simple form p\{v)p2{£t) ■ I n section 
12.31 we describe how operators with non-smooth symbols of this form oper- 
ate in semi-classical Sobolev spaces. Classically, norms of pseudodifferential 
operators depend on derivatives in x and in £ of the symbol, and, because 
x and £ play somehow symmetric roles, derivatives in £ can be shifted to 
derivatives in x. Here we can afford to lose derivatives in x. The eventual 
energy estimate in H£ involves ||ita||# s ' , with s' > s. This does not harm the 
proof if the initial datum is assumed to have enough Sobolev regularity. 

1.5 Background and references 

The (Z) system was derived from kinetic models by Vladimir Zakharov and 
his collaborators in the seventies jlSj . 

The initial value problem for the (Z) equations has received much at- 
tention. Global existence of smooth solutions in one space dimension (and 
of weak solutions in two and three space dimensions, for small initial data) 
was proved by Sulem and Sulem [221 ■ Global existence of smooth solu- 
tions in two space dimensions, for small initial data, was proved by Added 
and Added pQ. Schochet and Weinstein [21] and Ozawa and Tsutsumi |19j 
showed existence of local in time, smooth Sobolev solutions. Colliander and 
Bourgain [3] and Ginibre, Tsutsumi and Velo ^2j studied critical regularity 
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issues for local solutions. For large initial data, no evidence of singularity in 
finite time is known in space dimension greater than one. 

In their book on the Schrodinger equation |23| . Catherine Sulem and 
Pierre-Louis Sulem show how the Zakharov equation can be formally derived 
from the Euler-Maxwell equations; the WKB asymptotics of section IcS.ll is 
based on their description, and on discussions with Vladimir Tikhonchuk 
and Thierry Colin. 

To our knowledge, the results of and [E] were the first results estab- 
lishing rigorous links between Euler-Maxwell and Zakharov. 

Formal WKB expansions, carried out in [2S], have shown how the weakly 
nonlinear limit of (EM) fails to describe nonlinear interactions; such a phe- 
nomenon had been observed by Joly, Metivier and Rauch in the context of 
the Maxwell-Bloch equations. Joly, Metivier and Rauch's paper Jl], that 
describes large-amplitude solutions of semilinear systems of Maxwell-Bloch 
type by means of normal form reductions, is the main inspiration of the 
present work. 

In Klein- Gordon- waves systems were formally derived from Euler- 
Maxwell, and the Zakharov equations were rigorously derived as a high- 
frequency limit of these Klein- Gordon- waves systems. The stationnary phase 
arguments of [Bj, where solutions were represented, through Fourier analysis, 
in the form e _1 JjJ e***/ £ B(t')dt' , are analogous to the normal form reductions 
of the present work. The above integrals are bounded if the ratio B/Q is 
bounded, which echoes the transparency condition ()1.13j) . 

Highly-oscillating, large-amplitude solutions of quasilinear systems were 
considered by Serre in [2U] , and by Cheverry, Gues and Metivier in [S] . These 
papers deal with conservation laws, in particular, non-dispersive systems, 
unlike the Euler-Maxwell system. In 5 , Cheverry studies the parabolic 
relaxation of the instabilities put in evidence in jS] and applies his results 
to the equations of the large-scale motions in the atmosphere. 

In ^U], Grenier studies a class of singular equations of the form 
with A of the form (|1.3|) . and B = 0. He proves existence of solutions over 
time intervals independent of e, under the assumption that A possesses 
a 'good' symmetrizer, in the sense that no singular terms are created by 
subprincipal symbols occurring in the symmetrization process. Grenier is 
naturally led to study operators depending on x through v (x), where v has 
a Sobolev regularity. He does not assume that the initial data are well- 
prepared, and studies the high-frequency behaviour of the solutions. 

Lannes recently gave precise bounds for the norms of pseudodifferential 
operators depending on x through v(x), where v has a Sobolev regularity, 
and for the norms of commutators of such operators. These questions had 
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previously been considered by Taylor j^H] , and by Grenier in the article men- 
tioned above. We use a consequence of Lannes' description of the paradiffer- 
ential remainder (formulated as Proposition 12.11 it is used in sections 12.5.11 
and 12331) . 

In the approximate solution u £ a to (EM) that is constructed in section 
13.11 the envelope of the electric field is O(l), while the fluctuation of density 
has size 0(e), in particular, it vanishes in the high-frequency limit. How- 
ever, in the (Z) system, the fluctuation of density has a finite effect on the 
electric field. This means that there is a strong coupling between E and n 
in the Euler-Maxwell equations, an evidence of which is the non-transparent 
condition 1)2. 24|) . Such a phenomenon was called "ghost effect" by the Kyoto 
school of Sone, Aoki and Takata. These authors extensively studied this phe- 
nomenon in the context of small Knudsen number analysis of rarefied gases; 
a good reference is Sone's book j^J, and the paper [55]. As Sone explains in 
|21j . ghosts effects are characteristic of situations where large temperature 
variations are recorded. It would be interesting to show a formal similarity 
between their formal Hilbert expansions (describing the continuum limit) 
and the WKB expansions of highly nonlinear geometric optics (describing 
high-frequency limits). 

M. Colin an T. Colin propose a generalization of the Zakharov system 
in [7j. Their system consists in four Schrodinger equations coupled with 
quasilinear terms and a wave equation. It describes three- wave interactions, 
in particular, the generation of a Raman backscattered field. It is an inter- 
esting question, to know whether or not the result of this paper could be 
generalized to their extension of the (Z) system. 

We conclude this introduction by mentioning open questions and direc- 
tions for future work. 

It is natural to ask whether or not the result still holds when the initial 
condition is assumed to be oscillatory, that is, has the form a e (x)e tk / £ , 
where a £ is a bounded family in H s . It is shown in |28| that, if k ^ 0, the 
limit system is (Z) c , where c = uj'(k), and a; is a local parameterization on 
the characteristic variety (see figure |2J). Linares, Ponce and Saut prove in 
|16j that this system is well-posed in Sobolev spaces ^Hj; Colin and Metivier 
prove in [Hj that it is ill-posed in L°°. 

Another interesting direction for future work is to consider the case of 
large perturbations of the initial data, of the form e k °<f e , with ko < 3 + i. 
Our guess is that the strong coupling between the electric field and the mean 
mode of the fluctuation of density would then create strong instabilities in 
short time. 
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2 A class of singular equations 
2.1 Symbols 

We consider profiles u,v,. . . depending on x € M. d , with values in C n , and 
symbols p,q, . . . depending on {e, v, £) G (0, 1] x C" x R rf , or on e, x, £, with 
values in the n x n matrices with complex coefficients. 
For < e < 1 and s e R, we let 

\\v\\ £ , s := \\(l + H\ 2 y /2 mh*(niy 

In particular, || • ||i )S denotes the classical Sobolev norm. A profile v is said 
to belong to Hg(M. d ) when ||u|| es is finite. The space Hf(W i ), or simply 
H s (M. d ), is the classical Sobolev space. Remark that 

ll^lll.s = \\v\\e,s, 

where 

h E v(x) := e d/2 v{ex). 

For k G N, let 

\v\k,oo ■= ^2 SUP \ d " V \- 

o<H<fc* eRd 

A profile w is said to belong to W k '°°(M. d ) when | v\k oo is finite. Let do > 
For all keN, the embedding F e fc+ * ^ PF fc '°° holds: 

|w|*,oq < Ce~ fc - d/2 ||i'|| eife+do . 

We now define, and somehow adapt to our context, the class of symbols 
studied by Lannes in ^S] (see also Taylor \2jo\ and Grenier [TU]). 

A symbol p(s,v,£) defined in (0,1) x C n x R d , is said to belong to the 
class C^M 171 , m G N, when there exists £o > such that 

• P|{|*|<1} 6 G°°((0,eo) x m. d ,L°°({\Z\ < 1}), and 

• for all a, /3, there exists a non-decreasing function C a ^ such that for 
all v, 

sup sup (e) l/3| - m |a £ %afp(e, W ,0!<C ai/3 (|t;|). (2.1) 

£G(0,£o) |CI>V 4 
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In particular, if p G C°°M m , then for all a, d* v p G C°°M m . 

A symbol p G (7°°^™ is said to be k-regular at the origin when 

P|{|€|<1} e ^((O. £ o) x R d , W fc '°°({|£| < !})• 

A symbol is said to be smooth at the origin if it is /c-regular for all k. Symbols 
in C°°M m that depend only on £,£, are called Fourier multipliers. 

If p G C 00 ^™ is evaluated at v G s > \ , then Moser's inequality 
implies that 

SUp ||p(u(-),f) -p(0,Olle,s < C , a (|v|o,oo)||w||e 1 «, 

I5I<1 

and for all /?, 

SUp (1 + \ef hm )\\0^ (pKO.O -P(0,6)lle, s < C / J,.(Mo,oo)|M|.,., 
ICI>l/4 

for some nondecreasing functions C s ,Cp iS . 
2.2 Para-differential operators 

The class 5™, m G R, fc G N, is defined as the space of symbols q(e, x,£) 
such that there exists < £q < 1) such that, for all < e < £q, for all a, (3, 
with \a\ < k, there exists C £ , a ,/3(t) such that for all 

(^ a \- m \d«dlq{e,x^)\<C £ ^. 

With this definition, if p G C 00 ^™ is smooth at the origin, and if v is a 
profile in W k '°°(M. d )), then q := p(v) belongs to S™. 

To q G 5™, one associates the pseudo-differential operator op £ /(g), < 
sf < 1, formally defined by its action as 

(o Pe ,(g)*)(s) := (2^ d / 2 / e^g( £ , x, 

With this definition, 

op £ (g) := (K^op^he, 

where q(e,t,x,£,) := q(e, t, ex, £). 

Symbols in 5™ are smoothed into paradifferential symbols as follows. 
Let xo '■ R+ — *• R be a smooth function, such that < x < lj an d 

Xo(A) = l, for A < 1.1; xo(A) = 0, for A > 1.9. 
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For k > 0, define ip k : R d -> R by 

MO ■■= Xo(2- fc (l + |e| 2 ) 1/2 ) - Xo(2- (fc - 1} (l + |£| 2 ) 1/2 ). 
With these notations, for all £ 6 R d , 

fc>0 

Let x : R d -> R be a smooth function, < x ^ 1; an d such that 
x(iy) = 1, for \rj\ < 1.1; x( 7 ?) = 0> f° r I 7 ?! — 1-9- 
Let ^ : R d x R d -> R be defined by 



k>0 



Then ip satisfies, 



Hv,0 



1, |r ? |<2- 5 (l + |£| 2 ) 1 /2, 
0, \r]\ >2- 1 (l + |£| 2 ) 1 /2. 



In particular, 



^(»7,£) = 1, for alU, for all |r/| < 2" 5 . (2.2) 

One lets ^ 

g^(77,£) := ip(r},£)q(ri,€). 

The paradifferential operator associated with g is 

opj(g) :=op £ ,(^). 

The following proposition describes how well the action of a pseudo- 
differential operator is approximated by the action of its associated para- 
differential operator, a classical result in the case of differential symbols 
and when e = 1, of which Lannes gave an extension to pseudo-differential 
symbols in ^3]. We check below that the result of Lannes extends to < 
e < 1; in particular, that the action of the para-differential remainder, in 
is very small with respect to e, when s is large. 

C denotes non-decreasing functions, and do is a real number such that 

[f ] < do < [i] + 1. 
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Proposition 2.1 (remainder). Let v G H£, s > |, and p G C°°.M m 6e 
smooth at the origin. Then 

op e (p(«)) = op^(p( V )) + ^/ 2 op £ (^))> 
and, /or a// w G # e m+d °, 

\\op £ (R p (v))u\\e,s < C(Mo,oo) IMMMkm+do ■ ( 2 - 3 ) 



Above, C denotes a nondecreasing function, and do is a real number such 
that [£] < d < [f] + 1. 



Proof. We indicate how Q2.3|) follows from Propositions 3.3 and 3.4 of |15) . 
Let q := p(v) — p(0). The operation of para-differential smoothing is a con- 
volution in x, and thus R q = R p r v y The Fourier transform of the symbol of 
R q is 

R q = (l-i;(ri,0)q(v,0- 
The point is that because of (|2.2j) . the above symbol is identically zero for 
small n. Let </3o := X> an d for > 1, 

<p k (r,) :=x(2- fe r/)-x(2- (fc - 1) r 7 ). 

Then, for all 77 G M d , 

k>0 

One can write 

|fc-fc'|>3 |&-fc'|<3 

where the sums run over integers k, k 1 > 0. The first sum in q is further 
decomposed into qi + q2, where 

Qi'-=Y Y fik'ivMVjOVkit), & : =Y Y <Pk>(v)q('il,0<Pk(0- 

k>3k'<k-3 k'>3k<k>~3 
Remark that 

k<3 

Thus we have 

Rq = Q2 + q r ,l + qr,2, 
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where 

fc'>0,|fc-fc'|<3 

q r ,2 := x(r/)^,6xo(2- 3 (l + |^| 2 ) 1/2 )-^x(2- fc+3 r/)^,6^^ 

fc<3 

The symbols q2, Qr,i, Qr,2 correspond to the symbols an, 07^1 and aR,2 in 
[T5| . We want to bound 

||op e (i2 5 )u|| E , a = ||op a (g 2 + Qr,i + q r ,2)h s u\\i )S , 

where q(x) := q(ex). Propositions 3.3 and 3.4 of ^3] imply that 

llopife + q r ,i)h e u\\ 1>s < M(q)\\u\\ etm+do , (2.4) 

where 

M(q):= sup sup (1 + |£| 2 )M- m ||fl£?(-, Oik- 
M<70 £eR d 

where 70 depends only on d. Now, owing to (|2.2|1 . g can be replaced by 
(1 — op 1 (x))(Z in the symbol of R q , the function \ being smooth, identically 
equal to one for \rj\ < 2 -6 , and identically equal to zero for \rj\ > 2 -5 . For 
all w G .fP, for all |&| < s, there holds 

||(l-o Pl (x))^||i, s <C||^|| 1)S _ fc . (2.5) 

Applying (|2.5|) to (|2.4j) . one finds that the contribution of (72 + Qr,l to the 
operator norm of R q is bounded by 

Ce s - d ' 2 sup sup(l + |e| 2 ) l7| - m ||a^(pW-p(0))|| £ , s , (2.6) 

M<7o ?eR d 

which, in turn, is bounded by e s ~ d / 2 C{\v |o,oo)IMks- The support of q r ^. is 
included in a ball \rj\ + |£| < A. Thus we can use Lemma 4.4 of |15| to 
estimate the contribution of op £ (g rj 2). Using (j2.5|) again, one sees that, up 
to a multiplicative constant, it is also bounded by (j2.fi j) □ 

The next propositions are classical results of para-differential calculus 
(see for instance Appendix 2 of JZ1)> formulated in our C^Ai 171 framework. 
In the following statements, C denotes non-decreasing functions, and cLq is 
a real number such that [|] < do < [|] + 1. 
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Proposition 2.2 (action). Let v G L°° and p G C^M™ 1 be smooth at the 
origin. The operator opf(p(v)) maps H £ +m to H£, for all s, and for all 
u G H s £ +m , 

\\opf (p(v))u\\ E>s < C(|i;|o )0 o)||it|| e ,s+m- 

Proposition 2.3 (composition). Let pi G C°°Mf and p 2 G C°°M™ be 
smooth at the origin, and let 

p4pi := £ ^T-^Pl^Pa- 
a! ^ 

|a|=l 

If PI1P2 are evaluated at v±,V2 G VF 1 ' 00 , /or aZZ it G H £ +1 ~ mi ~ m2 , 

II (op^(pi)opf (p 2 ) - op^(pip 2 ))« || e , s < eC(|ui,w 2 |i,oo)|tt|| e , s +i —mi — m2 • 
If PI1P2 are evaluated at v\,v 2 G W 2 '°°, then for all u G H £ +2 ~ mi ~ m2 , 

II (op^(pi)op^(p 2 ) - Opf(pip 2 ~ epi$p 2 ))u \\ £ ,s 

< £ C(| V\ , V2 1 2, 00 J |^||e,s+2— mi— m2 • 

Proposition 2.4 (adjoint). Lei v G W 2, °° and p G C 00 ^" 1 be smooth at 
the origin. Let p(v)* denote the complex adjoint of the matrix p(v), and let 
°vt(p( v ))* denote the adjoint ofopf(p(v)) in L 2 . Let 

r *(p) ■= E ^r-^ P (vy. 

|o|=l 

Then, for all u G H s £ +m ~ 2 , 

|| (opi(p(v))* -opt(p(v)*) -eopt(rMv))))u \\ e>8 

< E 2 C(\v\2,oo)\\u\\ £ , s+m -2- 

2.3 Pseudo-differential operators with limited regularity 

We now consider non-smooth symbols that have the simple product struc- 
ture: 

p(e,v t £) =Pi(v)p 2 (e,£). (2.7) 

If p is a symbol in C°°.M m , with the structure (|2.7|) . then, in particular, p\ 
is a smooth map, and p 2 G C°°Ai m . Matrix- valued symbols are said to have 
the structure (|2.7|) when every entry can be written as a sum of terms of the 
form 1(277)1 . 

We will use the following lemma: 
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Lemma 2.5. Let u, v & H s (R d ), s > |, and assume that v has compact 
support. Then 

\\uv\\i,s < C\\u\\i tS \v\o j0 o, 
where C depends only on s and on the space dimension. 

Proof. This estimate follows easily from a dyadic decomposition. More de- 
tails can be found in Lemma 3.1 of |15j . for instance. □ 

In the following propositions, C denotes nondecreasing functions, and 
[#] < do < [|] + 1. 

Proposition 2.6. Letp G C°°M m , of the form (|2~7f) . andv G H s e , s > 
For all u G H s £ +m (W i ), 

||op e (p(w))w|| E)a < C(||«|| £)S+m + ||«||i,dol|«||i,2 + s s ~ d/2 \\v || e , s ||n|| ejTO+(io ), 

where C depends on |f|ioo- 

Proof. If p does not depend on v, the result obviously holds, and so, changing 
Pi to p\ — pi(0) if necessary, we are reduced to the case pi(0) = 0. 

We use the smooth truncation x introduced at the beginning of section 
12.21 As p2 is smooth for |£| > 1/4, the action of p\(v)(l — op e (x))op £ (p2) 
can be estimated with Propositions 12. ll and 12.21 It remains to bound 

\\0Y> £ {pi{v)xP2)u\\ e ,s = lbl(w)0Pi(XP2)fo e «||l, a , (2.8) 

where v{x) = v(ex). We can write pi(v) as the sum of op 1 (x)pi(u) and 
(1 — op 1 (x))pi(v), and apply the above Lemma. Thus we need to bound 

\op 1 (x)Pi(v)\o,oo\\op 1 (xP2)h £ u\\i, s , (2.9) 

and 

11(1 - oVi(x))Pi(v)\\i, s \ov(xP2)h E u\ 0iOO . (2.10) 
The first factor in (|2,9|) is bounded by 

sup | [ e^xieOn&mnl < C\\pMv)\\ L i 

X J 

< C(\v\ 0:OO )\\v\\ 1Ar 

With the first factor in (l2~TI71) is bounded by e s ~ d/2 C(|u| ,oo) ||u|| e , s . 

Finally, the bounds 

\\op 1 (xP2)h s u\\x tS < C||ti|| L 2, |op 1 (xP2)^ £ «|o,oo < C\\u\\ e4o , 

yield the desired estimate. □ 
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Next we describe the composition of two operators of the form (|2.7|) . 
Because the composition of two such operators involve remainder terms that 
do not have the form (|2.7[) . we need the following notations and lemmas. 

Let q be defined on R|, and p be smooth on M^. Introduce the notations, 

Q £ (£,0 + 0-9(0, 

and, formally, 

P(ff S P)|(e,*,0 := (0 £ (-,0 *p)( ic )» 

where Q £ denotes the inverse Fourier transform of Q e in its first variable, 
and * is a convolution in x £ 

With these notations, op e (/?(g,p))u is, formally, the inverse Fourier trans- 
form of 

J (q(sO - q(eO)P&)(£ ~ 0^'K'- 

Lemma 2.7. Let p be smooth on R n , such that p(0) = 0, and \ be the 
smooth truncation introduced at the beginning of section \2.2L 

(i) #9{|e|<l} andv^H s+d °, 

\\op £ {p{xq,P{v))u\\s,s < C(|t>|o,oo)||w||l, s+ d ||?x|| L 2 

(ii) I/g{|e|<i} G W 1 ' 00 ond « G #«+*>+!, then 

\\0Pe(p(xq,P(v))u\\ £:S < £(7(Ho,oo)IM|i, S +^+i|M|l2. 

(iii) If q £ C^M 171 is a Fourier multiplier, if v G jp+ao+m ^ ^ en 

ll°P £ (p((l ~ X)q,P(v)))u\\ £ ,s < eC(\v\ Q>00 )\\v\\ li8+m+do \\u\\ e)S+m -i. 

The third estimate is not tame, but it will be sufficient for our purposes. 

Proof, (i) Let go := xq, an d wq := op e (/9 £ (go ) p(^))'U- The fl^ norm of wq is 
equal to the L| norm of 

(i + H\ 2 yi 2 J (goH') - qoH))pV)^ - 0«(0^', (2.11) 

With Peetre's inequality and because go is compactly supported, the L| 
norm of (|2.11|) is bounded by 

Mo,oc ii j\i + w - oi 2 )' /2 ipSou - o\m'm'\\ LV 
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The above integral is a convolution. Thus 

IKIks < c\\t- 1 \{i + h\ 2 y' 2 ^v)\t- 1 \u\\\li 

< C\F-\l + \eey /2 l&)\\o,oc\\u\\ L 2. 

Now the L°° norm of JF-^l + \ei\ 2 ) sl2 \p{v)\ is bounded by the L 1 norm 

of (1 + H\ 2 ) s / 2 p(v), which in turn is bounded by C(|u|o,oo)||'y||i,s+(io- This 
yields the desried estimate. 

(ii) If d2qo is bounded, then w is the inverse Fourier transform of 

£ E / / d M £ t + - o) dt aSSte - em'n', 

ItI=i 

and the bounds that led to (i) are easily adapted to obtain (ii). 

(hi) Let q\ := (1 — x)<7, and w\ = op e (p((l — x)QiP( v ))) u - The H£ norm 
of w\ is equal to the L| norm of 

I7l=i J Jo 

There exists C such that, for all 

\dJqM + - 0)1 < + ke?) (m - 1)/2 + (i + - oi 2 ) (m ~ 1)/2 ). 

The rest of the proof of (hi) is similar to the proof of (i). 

□ 

Given two symbols of the form (|2.7|) : 

p(e, v, £) = pi(v)p 2 (e, £), q(e, v, f ) = q 1 {v)q 2 (e, £), (2.12) 

let 

*(p,q) ■= qip(xq2,Pi)p2- (2.13) 

Lemma 2.8. Letp,q of the form (j2,12|) . where pi,qi are smooth and vanish 
at v = 0, x?2 G and y>2 G C°°.M m . Let r := r(p,q). If p\ and q± are 
estimated at v E ^ or a ^ u ^ ^r, 
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||op £ (r)u|| e , s < C(||-u||i iS+ d )||n|| £im . 
and, if a £ C^M 1 is a Fourier multiplier, for all u S H™, 
ll(op e (o)op e (r) -op e (ar))u|| eiB < eC(||u|| M+(io+ i)|| 



u 
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Proof. Let z := op £ (p2)u, p\ = pi(v), and q\ = qi(v). The Fourier transform 
of op e (r)u is a convolution in £, 

m - e')((x«a)(^") - (xfe)(ec'))w(e' - o^")«, 

and, because XQ2 £ ^ 00 ; the first estimate in the lemma is obtained in the 
same way as Lemma 12.71 (i) . 

Let w := (op £ (a)op e (r) — op £ (ar))u. The Fourier transform of w is the 
sum, over I7I = 1, of 

e J a 7 ( £ , £ £, £ OS7l(£ " 0((X92)(eO " (x<? 2 )(^'))pi(£' " t'WW , 
and 

e / a 7 (£,^^ £ e / 0fl(e-0((x92)(ee / o-(xg2)(eeo)^^^(c , -n^(r)«^ 



where a 7 (e, 77,?/) := dja(e,i] + t(t]' — 7]))dt. Again, because a 1 ,XQ2 6 £ oc \ 
these convolutions can be bounded in the same way as in the proof of Lemma 
12 .7\ to yield the second estimate. □ 

We can now state a proposition that describes the composition of two 
symbols of the form (|2.7j) . 

Proposition 2.9. Given p G C°°X m i and g G C°°A4 m2 , 

(i) i/p and q have the form (|2.12j) . i/p is 1-regular at the origin, if p and 
q are estimated at v G # s +do+l ) / or a // u 6 ^s+m 1 +m 2 -l ) 

II [op e (p),op e (g)]n - op e ([p,g] +r(p,g))n || £)S 

< eC(||w||i iS+ d 0+ i)||n|| £iS+mi+m2 _i; 

(ii) if p is a Fourier multiplier and is 2-regular at the origin, if q has the 
form and is estimated at v G H s+do+2 , for all u G H'|+mi+m 2 -2 ) 

II [op e (p),op e (g)]n - op e ([p,g] + e(pjjg))u || £)S 

< e C(||'y||i )S+( i 0+ 2)||u|| ejS + mi + m2 _2- 

Proof, (i) We compute 



op £ (p)op e (g)n = Piop £ (p 2 )(qiop £ (p 2 )u) 

= Pi(op £ (P2Qi) ~ op £ (p(p2,qi)))op £ (p 2 )u. 
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Thus 

[ap e (p),op e (g)] -op e (\p,q]+T(p,q)) 

= Piop £ {p{p2,qi))op £ {q 2 ) -giop e (p((l -x)?2,Pl)op e (p 2 ), 

and the action of the right-hand side on it is estimated with Lemma l'2. 71 (ii) 
and (iii). 

(ii) Let w := ([op e (p), op e (q)} — op e [p, q] +eop e (p$q))u. The Fourier trans- 
form of w is 

where p 7 (ef,ef) := f (l-t)djp(e£' +et(£-£'))dt, and we obtain the second 
estimate as above. □ 

2.4 Assumptions and results 

We use the profile and symbol spaces introduced in section |2~T1 Let A, B , Q e 
such that 

• A is a smooth symbol in C°°ftA l ; 

• B is a bilinear map M n x M n -» M n ; 

• <7 £ is a family of smooth maps: M n — > R n , such that (/ e (0) = 0. 

Using the notations of section |2~TT we denote by op £ (A(e, ev)) the semi- 
classical pseudo-differential operator with symbol A(e,£V,£). We write the 
Taylor expansion of A in v as follows: 

A(e,ev) = A {0) (e, 0) + eA {1 \v) + e 2 A {2) (e,v), 

where A® G C 00 ^ 1 , for j = 0, 1,2, and where ^(u) := <%.4(0,0) • v is 
linear in u and has the form ()2.7jl . 
We study the initial value problem 



?/(0,x) = a £ (x) +e fc V £ (^), 



where the initial datum a e belongs to i? <T (R ), for some Sobolev index a, 

d 
2 



much larger than ^ (we will actually need a > 6 + d) , uniformly with respect 



to e : 

0<e<e 



sup \\a £ ||i o- < oo. 
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The perturbation e °(p £ is such that 

k >3 + ^, (2.15) 
and ip e belongs to H*(R d ), uniformly with respect to e : 

SUp ||^ £ || £ , s < OO, 
0<e<£ 

for some large Sobolev index s, smaller than a. 

Our first assumption is a hyperbolicity assumption that implies in par- 
ticular the local well-posedness of the initial value problem (|2.14j) . 

Assumption 2.10 (hyperbolicity). For alle,v,£, the matrix A(e , v , £) is 
hermitian. Let 

A= Yl iX J Il o+ Y iAfcIIfc ' ( 2 - 16 ) 

l£j< n no+l<fe<n 

be its spectral decomposition, where the eigenvalues \j,\k « r e real and the 
eigenprojectors IL;,IIfc are orthogonal. We assume: 

(i) the eigenvalues can be ordered as follows: for all e, v, £, 

sup \\k(e,v,£)\ < sup \\j(e,v,g)\; 

no+l<fc j<n>a 

(ii) for alll<m<n, X m (e, 0, £) G C^M 1 , and Yl m (s, 0, f ) € C^M ; 
(hi) for all uq + 1 < m < n, A m (0, 0, £) = 0, for all £. 

In Assumption 12.1(11 we do not assume that the eigenvector decomposi- 
tion is not singular, as we want to handle the case of the (EM) equations, 
whose eigenvector decomposition does become singular for small frequencies 
(see figure |U and section lS. 2.1(1 . 

Let 

n := Y n * := Yl U r ( 2 ' 1? ) 

l<fe<no no+l<j<n 

In addition to Assumption 12.101 we will assume that 

U ,U s (e,v,0 eC°°M°. (2.18) 

When A depends analytically on e, v, £, ()2.18|) follows from standard con- 
siderations, as detailed in section |HI 
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In reference to the Euler-Maxwell equations (see section |3J) , the eigenval- 
ues Xj for 1 < j < no are called Klein-Gordon modes, while the eigenvalues 
Afc, for no + 1 < k < n, are called acoustic modes. Condition (i) in As- 
sumption [2211 states that the acoustic modes do not cross the Klein-Gordon 
modes. Condition (ii) is a regularity assumption, and condition (iii) amounts 
to say that the acoustic velocities are 0(e), a consequence of (jl.7j) . 

Assumption 2.11 (approximate solution). For all l < a — 2 — |, there 
exists t*(lo) > 0, independent of e, and a family of profiles u £ a , such that 

dt< + ^op £ (^l(e, e<))< = <) + £ e «) + e l °R £ a , 

u £ a (0,x) = a e {x), 
where R £ a is uniformly bounded with respect to e, 

sup sup ||i^(i)||i )O ._z _ 2 < oo. 

0<e<e 0<t<t* (l ) 

There exists a characteristic frequency lo ^ 0, a finite set 1Z* C Z and profiles 
{ u a,p}p£Tl* an d "ai such that u £ a decomposes as 

u e a (t,x)= Yl e ip ^ £2 u £ atP (t,x)+ev £ a (t,x), 
pen* 

with the uniform bounds 

sup sup (\\u e (t)\\i jCr + \\d t u £ a p (t)\\i )a - 2 ) < oo, 

0<e<eo 0<t<t* (l ) 

sup sup ( IK (*) || + lk 2 ^o(*)lll,<r-l) < oo. 

0<e<eo 0<t<t* (l ) 

For the Euler-Maxwell system of equations, such an approximate so- 
lution is explicitly contructed in section EH under an assumption of well- 
preparedness for the initial datum. In the following assumption, 1Z* refers 
to the set of characteristic harmonics introduced in Assumption 12 . 1 fl 

Assumption 2.12 (resonances). There exists < q < c m < C m , such 
that the resonance equations in £ G R d and p,p' G TZ*, 

® jAp (e) := A i (e,0,0-A fc (e,0,e)+l>w = 0, (2.20) 
«W := A i (0,0,0-(P + P')w = 0, (2.21) 

are such that, 
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(0-0) 1 < j, k < no : i/ie solutions £,p of (|2.2Uf) e = 0, are located in 
the interval c m < |£| < C m ; outside this interval, $^^(0) is bounded 
away from 0, uniformly in £, for all p G 72-*; 

(0-s) 1 < j < no, no + 1 < & < n : /or e small enough, the solutions £,p of 
(|2.2U|) are located in the interval |£| < q; outside this interval, $j fe p (e) 
is bounded away from 0, uniformly in £ and e, /or aZZ p € 7£*; 

(0-0-s) 1 < j < no : /or a// € 7£*, VPjpp' is bounded away from 0, uni- 
/orn% in |£| G [0,c m ]. 

Next we state the assumptions that describe the interaction coefficients 
at the resonances. 

Introduce first the notations: 

B(u%) := B(ul-)+B(;ul)-A^(-)ul (2.22) 
V(u%) := (^'(O-^tO.OK,-), (2.23) 

where u £ a is the approximate solution given by Assumption 12.111 and, for 
all z, d v A^ 2 \0, u £ a ){u £ a , z) is the linear term in z in .4( 2 )(0,n^ + z)u £ a — 
.4> 2 )(0, u £ a )u £ a . The symbols B and T> depend on u £ a and of its e-derivatives; 
both belong to C°°M°. 

We assume that for some j < no < k, there exists p, £o and i] > 0, such 
that $j,fc, P (£o) = 0, and 

|n,(0,0)S« jP )n fc (0,0)| >n, uniformly in £ ~ (2-24) 

Inequality ()2.24|) means that the interaction coefficient HjBTLk is not trans- 
parent for resonances between Klein-Gordon and acoustic modes. 
Let 

P (e,u £ a ) -.= (n s ^ ))tino+(n s M( 1 ))no+n s ^ )(notjno)-(n s ^)no, (2.25) 

where the projectors Ho,II s are evaluated at (e, and is evaluated 
at u £ a . Let 

B r (e, <) := U S (B + d u U s ■ (e 2 d t u £ a ) + ep)U , (2.26) 

where p is evaluated at (e,u £ a ), B is evaluated at u £ a , and the projectors and 
their derivatives are evaluated at {e,£u £ a ). Remark that B T (e,0) = 0, and 

B r (e, <) = d u B r (e, 0) ■ < + £^73 r (e, 0) • «, u|) + 0(e 2 ), 

The linear term d u B T (e,0) ■ u £ a is the crucial interaction coefficient. The 
following transparency assumption states that it is sufficiently small at the 
resonances between Klein-Gordon and acoustic modes. 
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Assumption 2.13 (transparency). There exists Eq > and C > such 
that, for all < e < £q, for j < uq < fc, /or a// p € 72.*, 

|n fe (e,0)2?«)n 3 -(e,0)| <Ce, (2.27) 

and 

\H k (e,0)(d u B T (e,0) ■<)n j -( £ ,0)| < C(e 2 + |% )A>p (e)|), (2.28) 
uniformly in |£| < q, x € M d and t S [0, t* (Zo))- 
Introduce finally 

/ n Bu o \ / n ^i o 

^ • V o n s en s y ' m • V o n*A 

where A and the projectors are seen as symbols depending on (e,eu), and 
B as a symbol depending on (e, u). 

Assumption 2.14 (symmetrizability). There exists S, a smooth Fourier 
multiplier in C°°.M , such that 

-\\ u \\e,s < {°Pe(S)u,u) £ , s < j\\u\\ 2 s , 

7 

for all u G ff| and /or some 7 > 0, and 



-(5£ + (SET) e C°°>! , -^(iSA + (iSAY) e C°°M 



In the following theorem, u £ a is the approximate solution at order Iq, for 
some 3 + ko < lo < a — 2 — | , whose existence is guaranteed by Assumption 
12.111 t* = t*(/o) is its maximal existence time, independent of e, and s is a 
Sobolev index, such that l + ^<s<a — Iq — 2. 

Theorem 2.15. [/rider Assumptions [2~Tnl lOTT I2T2I l2~T3l and l2~Hl iaere 
exists a unique solution u £ G C°([0, to], H^(M. d )) to the initial value problem 
(|2.14|) . /or aZZ < to < t*\ there exists C > and eq > 0, swca that, for all 
< e < £ , for all < t < t*, 



sup \\(u £ -ul)(t)\\e,s<Ce k0 -\ (2.29) 

0<t<t 



In particular, 



m)\o,oo<Ce k °- 1 - d / 2 . 
o<t<t 



sup I (u £ — u~ 



In the error estimate (|2.29[) . C depends on a Sobolev norm of the initial 
data and on to. 
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2.5 Proof of theorem [2351 



In the proof below, we often drop the epsilons as we write u for u £ , u a for 
u e a , etc. We use the notations and results of section |2~T1 to describe symbols 
and operators. 

We start with (|2.14[) . Let Iq > ko + 3, and let u a be the approximate 
solution at order Iq given by Assumption 12.111 An existence time for u a is 
t* > 0, independent of e. 



2.5.1 The perturbation equations 

The exact solution u is sought as a perturbation of u a : 

u = u a + e k °- l u, (2.30) 
The symbol .A*- -* does not depend on v. This implies, with Proposition l2.il 

ap e (A(e,eu)) = opt(A(e,eu)) + e s+1 ~ d/2 op e (R AW(u) + eR A( 2) {£iU) ). 
The perturbation equations are 

dtu + ^opt(A(e, eu))u = ^op^(£)n + op £ (P)n + eR £ , 
u(0, x) = e(p £ (x). 
where B and T> are given by Q2.22JI and ()2.23() . and where 

R e :=e *o-*(B(u,u) + / (1 -t){g £ )"(u a + e k °- l tu) ■ {u,u)dt) 
Jo 

-e^- 1 [ {l-t)d 2 v A {2) {u a + £ k tu) ■ (u a ,u,u)dt 
Jo 

-e^-^op^R^^ +eR A(2){£>u) ) -e l °- k °- 2 Rl 

Under Assumption 12.101 standard hyperbolic theory provides the existence 
of a unique solution u to (|2.31D over a small time interval [0, i*(e)]> with the 
uniform estimate 

sup sup ||w(t)|| E)a < 5. (2.32) 

0<e<eo 0<t<t»(e) 

The term R e is a remainder, in the sense that its i?| norm can be 
bounded in terms of 5, uniformly in e. The H£ norm of the terms in the first 
line in the definition of R £ is indeed bounded by e k °~ 2 ~ d ^ 2 C\\u\\ et d \\u\\ s ,s', 
the terms in the second line R £ are bounded by 

e ko ~ 1 ~ d C(\\u a \\ £ , s+1 \\u\\ 2 tdo + |« a |l,(»N|e,<iolH| £ , s ), 
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and the last terms are bounded by 

^- 1 - d / 2 C||n|| £ , 1+(io || U || e , s +^- fe »- 2 ||^|| £ , s . 

In all these estimates, C depends on |ti|o,oo> 

With the estimates for u a given in Assumption 12.111 the bound for ii 
given in (|2.32|) . and the form of the equation (|2.31l) . for a + |/3| < 2, 

sup \(e 2 d t ) a dP U |o,oo < c a + e k °- 3 - d / 2 5, (2.33) 

0<t<t t (e) 

where c a does not depend on e. The size of the perturbation of the initial 
data in (|2.14|) . namely 0(e k °) in H^, where ko satisfies (j2.15|) . was chosen 
in order that the estimate (|2.33f) be uniform in e. 

With these notations, the above estimates give, for £q small enough, and 
s> 1 + 

\\R £ \ks<CMs,s + e l0 - k °- 2 \\R £ a \\ £yS , 

where C is a nondecreasing function of £q° 3 d ^ 2 5, ||n a ||e,s+i; s and d. 

We generically denote by i?( ) any pseudo- or para-differential operator, 
possibly depending on the solution u, such that, for all z € (M d ), 

\\R{o)z\\ e , s < C\\z\\ e>s , (2.34) 

uniformly in t € [0,t*(e)], where C is nondecreasing, and 

C = C{e k Q ~ Z ~ d/2 5, ||n o || 1(S+d0+2 ,s,d), (2.35) 

and where [|] < d < [|] + 1. 

We denote by 0(e k ) symbols associated with pseudo- or para-differential 
operators of the form e k R^ y 

In the next section, we are led to study compositions of para-differential 
symbols of the form Pj(u a ^), where pj G C^M^ is smooth, and u a ^ = 
{e 2 dt) a dxU. It follows from Propositions 12.21 and 12.31 that 

• If m,j < 0, and a + < 2, then opt(pj (u a ^)) = -R(o)- 

• If mi + m 2 < 1, if pi = pi(u a ^), and p 2 = P2(u a ',P')> witn a + a ' + 
\P\ + W\ < 1, then 

opt(pi)opf (p 2 ) - op^(pip 2 ) = e-R(o), 

• if m\ + m 2 < 2, if p\ = p\{u) and p 2 = p 2 (u), then 

opt (pi)opf (p 2 ) - op^(p lP2 ) - eop^iflpa) = £ 2 #(o)- 
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2.5.2 Projection and rescaling 

Let Xe '■ — * ^ be a smooth function such that e < Xe < 1> an d such that 
Xe is identically equal to 1 for |£| < co, and identically equal to e for |£| > ci, 
where < co < q < ci < c m , with the notations of Assumption 12,121 
Introduce the change of variables 

v : = op^(il )n, v s := ^opf (x e )opf(!I s )u, v := (v ,v s ), 

where IIo,n s are evaluated at (e,e(u a +£ u)). Then 

u = v G + eo V t{Xe 1 )v s . (2.36) 

With Proposition 12.31 and the orthogonality of LTo and H s , 

opf{n )v = v + e 2 opf{n mo)u + e 3 R {0) u, 

opt(U s )v s = v s + eopf (x £ (n s Jjn s ) + (x e tJn s )n s )n + s 2 R {0) u. 

We multiply (l2~3H) by op^(n ) (resp. op^(x e )op^(n s )) to the left to find 
the equation satisfied by vo (resp. v s ). We use Proposition 12.31 to spell out 
the compositions. 
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The terms in dt are 



qp?(no)$« = dm - o P f(d t n ; 



u 



and 



d t v - -opf((edtILo)Uo)v 
-op^(ex~ 1 (e5 t n )n s )v s + J R( 0) 'y; 



^opf(x £ )opf(n s )9in = a tUs - ~opt(Xs)opt(d t Il s )u 

= d t v s - ^opt((ed t n s )n s )v s 
1 



^optUe(ed t U s )U )v 
- ^opt((xe^d t U s )U )v + R (0) v. 
The terms of order one are, with Assumption 12. 1(JI 
op^(n )op^(^)n = (opf(U A) + e 2 op^(n M) + e 3 R {0) ) 



u 

op^(n ^)i;o + eoptHex^f^AixMs)^ 

and 
1 



£ Opt( X e)opt(U s )opt(A)u = opt(IL s A)v s + eopt(ex7Hxd(nsA)))v s 

t(p s0 )vo + eoptl 



+ eopt( Ps0 )v + eopt(xe^sA)U )v + 



+ E 2 R(p)V\ 

where 

Ps0 ■.= Xe((n s A)mo + (n s M (1) )no + n s ^(n tjn ). 

In the above symbolic computations, the projectors Ilo and H s , as well as 
A, are evaluated at (e, eu a + e ko ii), and is evaluated at (e, u a + s k °~ l u). 
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The singular terms in the right-hand sides are 

o V t(U )op £ (B)u = optiUoBUo^o + opfiex^IloBU^Vs 

-opt(xe)opt(U s )opt(B)u = -opt( X£ U s BU )v 

e e 

+ o P ^(x £ (n s ^)no + x4(n^)n )t;o 

+ eR(Q)V. 

In the above symbolic computations involving B, the projectors ilo,n s are 
evaluated at (e,eu a ), and B is evaluated at u a . 
The other singular source term is 

^opt( X e)opt (U s )opf (V)u = iop^x £ n s OTL>o + R(p)V, 
where the projectors are evaluated at (0,0), and T> is evaluated at u a . 
The equation in v thus takes the form 

d t v + ^opf(iA)v = ^opf(B)v + ±opf (D)v + R (0) v + r £ a , 

where 

• A = A(e,e(u a + e^" 1 ^)) e C^M 1 is defined as 



iA := 



AU 

o au s 



B = B(e, u a ) e C°°M° is defined as 



B :-- 





R sQ 



where 



B s0 :=Xe(n s (B + eV)U + e(d t U s )U + e(p s0 - (n,^)n ) 
-exe^ed t U s + U s A-U s BU ). 

where IIo,n s ,^l are evaluated at (e,eu a ); A^ is evaluated at (e,u a ), 
and B and V are evaluated at u a . 
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D = D(s, u a )=g + &' + E + F£ C°°M is defined as 
r> ._ ( M r" ._ ( £0, 

- ' V S* J ' - • V o o 

where 

:= - £X r 1 (x e tt(n s (^--B))n S) 

#0* := % H^oBUs + e9 t n - (ex^^^IT))!!, 



E := 



n £n 
o itsil 



where Ilo,n s ,^4 are evaluated at (0,0), /3, £> are evaluated at u a , and 

/ (^n )n \ 

- ' V {ed t U s )U s J ' 

where edtHj is short for <%IIj(0, 0)-(e 2 dtu a ), j = 0,s, and where Ilo,n s 
are evaluated at (0,0); 

. rl := e l °- k °-\opt(IL )Rl ±op^(IT)^). 
Next we polarize the source terms, by letting 

B s o ■= n s i? s0 no, b s := il s b s ii s , b 0s := n 0; B 0s n s , 

and 

"-(LI)- B "-(o B o* 



We let also 



and 



F + ( ^os - B o 
D:= B' + B" + E + F. 



The equation is now 

d t v + ^op^(iA)v = ^opf(B)v + iqp#(Z>)v + R {0) v + r%. (2.37) 

In l|2.37j) . the variables vq and v s are coupled only by order-zero terms, the 
leading singular term is polarized and has a nilpotent structure. The system 
is prepared. All the symbols in (|2,37|) are smooth. 

In the next sections, the terms B" , B' and B will be eliminated by normal 
form reductions. In the subsequent H£ energy estimate, the term E will be 
symmetrized, while the non-polarized term Fje will be seen to contribute 
toO(l). 
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2.5.3 First reduction 

In this section, the non-resonant term B" is eliminated from (|2.37|) . 

Proposition 2.16. Under Assumvtions \'2.V2\ and XL. 131 there exists a smooth 
symbol L G C°°M~ l such that 

[e 2 d t + opf (i A), opf (L(u a )] = opt(B")+eR {0) . (2.38) 

Proof. The leading term in the symbol of the source term in (|2,38|) is linear 
in u a : 

B 0s = d u B 0s (0,0) -Ua + Oie), 
and satisfies n 5 0s II s = B 0s . We look a solution L G C^M^ 1 to (j2~35|) in 
the same form: L = f ^ , where HqLq s II s = Lq s , and 

pen* 

where Lq SjP is linear in u ayP . Then 

[opt A), opt (L)] = [opt(iA(Q,Q),opt(L)}, 
up to the commutator 

[opt(A-A(0,0)\opt(L)], 

which, because A depends on u through eu, and because L is assumed to 
belong to C°°A4 _1 , has the form sR(oy If we suppose in addition that L is 
smooth at the origin, we can use Proposition 12.31 to obtain 

[opt(iA(p,0),opt(L)] =opt[iA(0,0),L] +sR {0) . 

Thus to solve (|2.38|) . it suffices to solve the equation 

e 2 d t L 0s + [iA(0,0),L] Qa = d u B 0s {0, 0) • u a , (2.39) 

up to 0(e). We compute 

e 2 d t L 0s + [iA(0 1 0) 1 L} 0s = Y, Y e* wt/ca ^,fc 1 p(0)n i L . iP n fc , 

j<no<k p&Z* 

where the projectors Uj,Hk are evaluated at (0,0), and where ^j t k, p (0) 
stands the evaluation at e = of the phase defined in (|2.20|) . Let xl be a 
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smooth function on M. d , identically equal to for |£| < q, and identically 
equal to 1 for |£| > cq. Let then 

L 0s , p := XL Yl ®ll P (°) U j( d u B 0s( > °) • ^, P )n fc . 

j<no<k 

Because &jl p i f° r 3 ^ n o < is uniformly bounded for j£| > q, the above 
defines a symbol L £ C7°°A^ — . Besides, this symbol is smooth, and solves 
(|2.39|) . up to the error term (1 — Xl)Bq s . Because \£Xs 1 Xl\ < £j this error 
is 0(e). 

□ 

Proposition 12.21 allows to evaluate the action of opf(L). Its norm, as 
an operator from i?| to is bounded by C(|u a |o,oo)- Consider now the 

change of variables 

v := (Id+eop^L))- 1 ?;, (2.40) 

Then 

8 t v = (Id+ eopt(L))d t v + opf (e^L)6. 
Because L is order —1 and is smooth at the origin, 

(Id+eopf(L))- 1 opf(A)(Id+eop^(L)) = op* (A) + e[op*(A), op^(L)] 

= opf(A)+£op*L4,£], 
up to error terms of the form e 2 i?( )- Similarly, 

(Id+eop^(L))- 1 op^( J B)(Id+eop^(L)) = opf(B) + eopf[B, L] + e 2 R {0) , 
and 

(Id+eop^(L))- 1 op^( J D)(Id+eop^(L)) = op? (£>) + e22 (0) . 
The leading term in e in the commutator 



[B,L] 



-LB s0 
B s0 L 



is 0(x £ )0(ex £ ) = 0(e). Thus the equation satisfied by v is 



$5 + ^opf(iA)v = ^opf(B)v + -opf(D)v + R (0) v + (2.41) 



i)v = -^op 
where f„ := (Id+ eopt {L))" 1 r £ a , and 



f»:=£)-[e 2 5 t + a-B,L], 
With the above proposition, 

D = B' + E + F. 
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2.5.4 Second reduction 

In this section, the non-resonant term B' is eliminated from (|2.41|) . 

Proposition 2.17. Under Assumptions \2A2\ a,nd \2.\?\ there exists a smooth 
symbol M G C^M" 1 such that 

[e 2 d t + opf (i A), opf (M(u a ))] = opt(B')+eR {0) . (2.42) 

Proof. The leading term in the symbol of the source term is linear in u a : 

B s = d u B s (0,0)-u a + O(e), 



and satisfies H S B S U S = B s . We look for a solution of (|2.42|) in the form 



M s 



of a smooth symbol M £ C°°AA. 1 , such that M = ( I , where 



n s M s n s = M s , and 

M S =Y^ e (ip "' )/£2 M s , p (n a , p ), 

where M SjP is linear in u atP . We check as in the proof of Proposition I2.1b1 
that in order to solve (|2.42|) . it suffices to solve up to 0(e) the equation 

e 2 d t M s + [iA(0,0),M} s = d u B s (0,0) ■ u a . 

Let 

M StP := ^ $7^(0)^(^^(0,0)-^)^, 

no<j,k 

where the projectors are evaluated at (0,0). Condition (hi) in Assumption 
I2.10l implies that for j, k > no, &jl p (0) = ipuj. The support of B s is included 
[co,ci]. Thus M SiP is a smooth symbol with compact support, M solves 

Proposition 12.21 allows to evaluate the action of opt(M). Its norm, as 
an operator from to H^ +1 , is bounded by C(|w a |o,oo)- Consider now the 
change of variables 

v := (Id+ eopt{M))~ l v, (2.43) 
With the above proposition, the equation satisfied by v is 

d t v + ^opf(iA)v = ^opt(B)v + ^opt(E + F)v + R (0) v + f £ a , (2.44) 

where f e a := (Id+ eopt (M)) _1 f^, and 

5 :=B + e[B,Af|. 
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2.5.5 Third reduction 

In this section, the resonant term B is eliminated from 1)2. 44 j) . under the 
transparency condition (|2.28|) . 

Proposition 2.18. Under Assumvtions 12.121 and |2~T51 there exists N € 
C°°M~ l , such that 

[e 2 d t + opt(iA)-eopt(E),opt(N(u a ))]=opt(B) + e 2 R (0) , (2.45) 



Proof. The source term is B = ( ~? ) , where 



B s0 

B s0 = Xe(B r + en s Pn ) + e(B™ - XeM s B T ) + 0(e 2 ). 
The symbol B r is introduced in (|2.26[) . and 

B^(e, Ua ) ■.= -n s ( X£ $(ed t n s + u s A w - n s £n ))n . 

With the transparency assumption (|2.27|) . the term ex e Tl s T)Ilo is 0{e 2 ). The 
Taylor expansions of B r and B nr in their second variables are 

£ r (e, u a ) = d u B T (e, 0)-u a + sd 2 u B*(s, 0) • (u a , u a ) + 0(e 2 ), 

and 

B* T (e, u a ) = d u B™{e, 0) • u a + 0(e). 
Thus, up to 0(e 2 ), the source i? s o is the sum of a linear term in u a , 

Xed u B T {e, 0) • u a + ed UJ B nr (e, 0) • u , (2.46) 

and of a bilinear term in u a , 

8Xs(d 2 u B T (0, 0) • (u a , u a ) - M s (u a )d u B T (0, 0) • u a ). (2.47) 

All the terms in (|2.46j) - (|2.47|) have the product structure (|2.7|) . Accordingly, 
we look for N in the form, 

N = iV (0) + e(N w +N {2) + A (3) ), (2.48) 

with 

JV(J) = [ J } ] G C 00 ^" 1 , for all i, 

and where 

42 



(a) iV (0) , iV (1) and have the structure l(2~7j) : 

(b) all the entries of have the form a(£)r(p, <?)&(£), for some Fourier 
multiplier a and b, and some symbols p, q in the form (|2.7|) (the nota- 
tion r is introduced in Q2.13|) ): 

(c) iV(°) and iV« are linear in u n , and and are bilinear in u n : 

N (J) = £ e^ 2 N^(u a , p ), j = 0,l, 

We now describe the symbols of the commutators 

[opt(iA-E),op E (N^)}, j = 0,l,2,3, 

using the results of section 12.31 

Because the symbol is assumed to belong to C^Ad^ 1 and to depend 
on x only through u a , 

[o V t{iA),o Ve {N^)\ = [op £ (iA(e,eu a )),op £ (N^)]+e 2 R m , 

and the commutator in the right-hand side of the above equation is 

[o Pe (L4M)),op E (iV( ))] +e[op £ (id v A(0,0) ■ Ua ),op £ (N^)\, 

up to £ 2 i?( )- Because N^°> is assumed to have the form (j2.7|) and A(e,0) is 
smooth, Proposition 12,91 (ii) implies that 

[op £ (iA(e, 0), o Pe (iV( ))] = op E [iA(e, 0), iV(°)] + eop £ (iA(0, 0)$N^), 

up to a remainder in e 2 . Because depends only on u a , Proposition 12.91 
(ii) implies that this remainder has the form e 2 R(Qy Remark that 

E{s,u a ) = d u E(0,0) -u a + 0{e). 

The symbol id u A — E is smooth, and thus Proposition 12.91 (i) implies that 
the commutator 

[op £ (d u U(0,0) • u a - E),op e (N^)}, 

is equal to 

opMiA " ^)(0,0) • u a ,N<®] +o Ve r{id v {A - £)(0,0) • u a ,N^), 
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up to a remainder in O(e). With Proposition 12.91 (i), and because iV"(°) 
depends only on u a , this remainder has the form e 2 R^y 
Similarly, for j = 1,2,3 

[op e (iA),op e (N®)] = [op e (a(0,0)),op £ (iV0-))] +eR i0) . 

For j = 1,2, Proposition 12,91 (ii) implies that 

[op £ (M(0,0)),op £ (ivW)] = op £ [iA(0,0),N^\ +sR (oy 

Lemma 12.81 implies that 

[op e (M(0,0)),op E (A^)] = o Ve [iA(0,0),N^} + eR {0) . 

These symbolic computations show that in order to solve (|2,45|) . it is 
sufficient to solve the system, 

e 2 d t N^ + [iA(e,0),N^U = Xed u B T {e^)-u a , (2.49) 

£ 2 dtN$ + [iA(0,0),N(% = d u B™(0,0)-u a - ^(0,0)^,(2.50) 

e 2 d t N$ + [iA(0,0),N^] s0 = B hl (u a ,u a ), (2.51) 

e 2 d t N$ + [iA(0,0),N^] s0 = r(d v (iA-E){0,0)-u a ,N^) s0 ,(2.52) 

with the notation, 

B hl (u a ,u a ) :=x £ (^ r (0,0) • (u a ,u a ) - M s (u a )d u B r (0,0) ■ u a ) 
+ [d v (iA-E)(0,0)-u a ,N(% - 
We now solve (|2.49jl . For j < n < k, for < e < eo, let 

*3,kA e ) - \ e 2/ 2) otherwise. 

These new phases are not continuous, but they are bounded in £, uniformly 
in e, for |£| < q. Let then 

N sZ-=Xe Yl ^j,kA^' lll k(e,OKd u B T (e,0)-u a:P )U 3 (e,0). (2.53) 

j<n <k 

Then every entry of is a product p(u a )q(e,^), where p is smooth, and 
the transparency condition (|2.28j) in Assumption 12.131 ensures that q(e, £) 
is bounded in £, uniformly in e. Besides, by definition, q is compactly sup- 
ported. Thus 7V~(°) G C°°M~ l , and it solves (l2~m 
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In (|2.50l) the source term A s (0, 0)^N^' is identically zero (Assumption 
EH(iii)). Let 

n $p : = E QjAP&y^kioMduB^M - Uai p)n 3 -(o,o), 

j'<no<fc 

Because S nr is supported in |£| € [co, ci], far from the Klein-Gordon/acoustic 
resonances, the above defines a symbol in C°°.M _1 . Then has the 

structure Q2.7|) . is linear in u a , and solves (|2.50|) . 

Let xn be a smooth truncation function on R d , identically equal to 1 for 
|£| < ci, and identically equal to for |£| > c m . Let then 

n $ p ,p> --=xn Y, *i y n s (o,o)i? bl K, p ,n ay )n J (o,o). 

The phases ^kpp' are uniformly bounded away from for |£| < c m (As- 
sumption ^. 121 (0-0-8)), Thus the above defines a symbol in C°°.M _1 . Then 
has the structure (|2.7|) . just like the source term B hl , is bilinear in u a , 
and solves (|2.51|) . up to the error (1 — XN)B bl . Because |(1 — Xl)Xs\ < e, 
the error is O(e). 
Finally, let 

^ioL' E ^^ y n s (o,o)rK iP ,n ay )n J (o,o), 

j<no 

with the notation, 

*{ua,p,Ua, P >) := r(d v (iA - E)(0,0) ■ u a , p , A^ 0) (n a y)) s0 - 

All the entries of 7V"( 3 ) have the form a(£)r(p, q)b(£), where a is smooth and 
p, q have the form (|2.7|) . The symbol iV"( 3 ) solves (|2.52)1 . up to the error 
(1 — XAO r > which, because is 0(e) for |£| > c\, is 0(e) as well. 

Finally ./V defined by (|2.48j) satisfies the assumptions (a), (b), (c), on 
which the symbolic computations were based, and solves (|2.45|) . 

□ 

Proposition 12,61 and Lemma 12.81 imply that the norm of op e (iV), as an 
operator from H£ to i?| , is bounded by O(||ii a ||i jS+ d 0+ 2). Consider now 
the change of variables 

w := (Id+ op e (iY)) _1 u, (2.54) 
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Then 

dtv = (Id+ eop £ (N))dtw + op £ (edtN)w. 
Because N is block triangular, 

(Id+ op £ (iV))- 1 op^(A)(Id+ op £ (iV)) = opf (A) + [qp^),op e (J\r)]; 

(Id+ op e (iV))- 1 opf( J B)(Id+ op £ (iV)) = opf(B), 

and 

op e (iV)op E (£d t iV) = 0, op £ (N)opf(E)op £ (N) = 0. 
Thus, with the above proposition, the equation satisfied by w is 



d t w + ^opf (iA)w = -opf(E + F)w + R (0) w + r £ a , (2.55) 

where r £ a : = (Id+ op e (N))~ 1 r £ , and 

opf(F) := (Id+op £ (iV))- 1 op^(F)(Id+op £ (iV)). 

2.5.6 Uniform Sobolev estimates 

We perform energy estimates on (j2.55j) . using the symmetrizer S whose 
existence is granted by Assumption 12 .141 
We evaluate 



d t (op £ (S)w, w) EiS = (op £ (S)d t w,w)e tS + (op £ (S)w,d t w) £ 
= 2M(op £ (S)d t w,w) £ , s . 

We can apply Proposition 12.41 to find, 

opt(iSAy = opt((iSAy) + e 2 R {0) , 
because A depends on u through eu. Similarly, 

opt(iSE)* = opt((iSE)*)+eR i0) , 

and this implies, 
1 

£ 

-M(op £ (S)opf(E)w,w) e , s = (op £ (-(iSE+(iSE)*)w,w) 

+ e 2 (R( )W,w) £:S . 



2®(op e (S)opt(iA)w,w) e , a = (opt(-(iSA + (iSA)*)w,w) £!S 
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Now Assumption 12.141 implies that the symbols in the right-hand sides of 
the above equations all belong to C°°M , and Proposition 12.21 implies that 
these symbols have the form R<oy The other source term contributes to 

^(op £ (S)opt(F)w,w) £ , s . (2.56) 

Introduce the notation, 

U:=( E ° ° V 
It follows from the definitions of the above changes of variables that 

w = op^(II)u> + eR(q)W, 
and that, up to a term of the form sR(p\, 

o P f(n)(id-op £ (iv)) = (id-o Pe (iV))opf(n). 

Besides, up to a term of the form eRr \, 

o V f(U)op £ (S) = op £ (5)opf (H). 
Thus, up to a term of the form e(Rr )W, w) E>s , (|2.56|) is equal to 
R(ap e (5)(Id- op £ (iV))op £ A (F)(Id- op e (N))w,w) e8 , 
where F := ILFII. Because IIo and IT are projectors, 

n ^n n = o, u s d t u s n s = o, 

and it follows from the definition of F (given in section l2,5.2j) . that F = 0(e). 
Gathering the above estimates, we find that 

d t (o P f(S)w, w) £ , s = (R (0) w, w) EiS + 2M(opt(S)f £ a , w) EjS . (2.57) 

The changes of variables of sections 12.5.31 to 12.5.51 define a normal form 
ty e (u a ), such that, 

w = (Id + ^ £ (u a ))v, 

and, for eo small enough and < e < eq, both ld+^ e (u a ) and (ld+^f £ (u a )) 
are uniformly bounded as operators H £ — * H £ , with norms depending on 
||w a ||i s+do+2- From (|2.57jl . we obtain, 

\Ht)\\l s <cy% s + c A||«(OII«,. + e , °-* - 8 l|ie(*')MIKt')lle,.cft', 

J 
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where the constant C depends on e ° 5, \\u a \\i tS +d +2, s,d, and on 7. 

Because (|2.36|) implies that u = R^v, Gronwall's lemma finally yields, 

\\u(t)\\ E , s < dW^lU^. (2.58) 

A classical continuation argument shows that, for sq small enough, the 
bound (|2.58|) is valid over a time interval [0,to]; independent of e. Then 
(|2.30[) yields the asymptotic estimate (|2.29|) . 



3 Application to the Euler-Maxwell equations 

We show in this section that the Euler-Maxwell equations satisfy the as- 
sumptions of Theorem 12. 151 The system we consider is (EM), introduced in 
section [T~2l in the specific regime Ql. 7|l - l|l. 8|) . For the unknown 



n; 



u £ (t,x) = (B,E,v e ,n e ,Vi, — ), 

a 



the system takes the form 



d t u £ + -^A(s,eu £ ,£d x )u £ = -B(u £ ,u £ ) +g £ (u £ 



where 



A(e,eu,ed x ) = Ao(e,ed x ) + eAi(e,u,sd x ), 
with the notations, 



(3-1) 



(3.2) 



A (e,0 :-- 



( 
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-ex 





i 















—i 





































i— 











eat; 




V 
















) 


diag(0 C 3 






■0,t 




e(vi 


■ I 



and 



B(u, u) := (0 C 3 , n e v' e , -9 e v' e x B, C , C 3 , C ), 
£ (u) := ^(0,f(n e )v e - ^{an i + ef £ {an i ))v i ,^^v i x B,0), 



with the notation, f £ (x) := e 2 (e £X — 1 — ex). 
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3.1 WKB approximate solution 

Let the initial datum 

a = (0,E°,u°, 0,0,0) G H a , (3.3) 
for some large a, where the electric field satisfies 

V • E° = 0, 

in accordance with 1|1.4|) . We assume that a is polarized (or well-prepared), 
in the sense that 

E° = E+(EY, v° = -E--E*, 

U> UJ 

for some fundamental frequency uj, defined in terms of uj pe , and some com- 
plex amplitude E (above, (E)* denotes the complex conjugate of E). 

Because the conservative form of the convective terms in the equations 
of conservations of mass in (EM)" allows simple formal computations, we 
carry out the WKB expansion on (EM)" rather than on (EM). 

Consider (EM)" in the regime 1)1.7(1 . We look for an approximate solution 
u £ app in the form of a profile 

u e app (t,x) = e[u £ app (et,x,9)] e=ojt/£ , (3.4) 

where u £ app has a WKB expansion, 

M 

u e \ e m u £ 

o-PP / j mi 
m=0 

such that for all m, u m is a trigonometric polynomial in 9, 
u m (t,x,0)= ^2 e ipe u m , p {t,x). 

The sets lZ m C Z arc finite <mcl th.6 v>m,p scr^ assumed to satisfy 

for some i*(e) > and some large k(m). We plug this ansatz in (EM)" and 
find a cascade of WKB equations, which we now describe. We sometimes use 
below the notation (v) p to denote the harmonic in 6 of a trigonometric 
polynomial v(i, x, 9). 
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Equations for the terms in 0{l/e 2 ). 

ujdgEo = v e0 , 
ivd e ve0 = -E , 
^(B ,n» ,v,o,n» ) = 0. 
The dispersion relation is 

u? - 1 = 0. 

We choose oj = 1. With this choice, and because e was set to be equal to 
(tOpeto) -1 , we find that e %wt l e = e iuipeT , where T represents the physical time, 
t = et()T, and oj pe is the electronic plasma frequency (|1.6|) . Thus the waves 
we consider are oscillating at the electronic plasma frequency. These waves 
are called plasma waves in the physical literature. The set of characteristic 
harmonics is 

K* = {-1,1}. 
The first term of the expansion satisfies 

£o,o = 0, v eQfi = 0, B = -B 0i0 , n^ = n S e0fi , v i0 = v m , nf = nf 00 , 

and 

E = E 0>1 e i6 + E ^ie- ie , v e0 = -E 0)l e ie - -E^ x e~ ie . (3.5) 

UJ UJ 

Equations for the terms in 0(l/e). 
udeBi + V x E = 0, 

a 1 

wd e Ei - V x B = v e i + n|l v e o - T-v i0 , 

ujd e v el = -6» e Vn^ - (Ei + 6 e v e0 x B ), 
uden^ + 6> e V • v e0 = 0, 

udevn = tJ-E , 



Let 



Lodgll^ = 0. 
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The frequency uj was chosen so that det (ipu + Lq) = 0. In C 6 , 

{ipuj + L )a = f ^ 

implies the compatibility condition 

ipuubi + &2 = 0. 

The oscillating terms satisfy 
and this implies, 

ipuj(nl v e0 ) p - # e (v e0 x B ) p = 0. (3.6) 
Besides, we find the polarization conditions, 



ipuB ltP = -V x E 0tP , ipujv iltp = —E , p , n\^ p = 0, n\ lp = j—^ V • E 0p . 
The nonoscillating terms satisfy, 



1 a 

- V x B ,o = v e ifi - —v i0fi: E lfi = e Vnl Q . (3.7) 

Equations for the terms in 0(1). 
iod e B 2 + d t B + V x Ei = 0, 

ujd e B 2 + dtE - V x Bi = v e2 + n^v e o + n|[ v e i - ^-(v a + nLv i0 ), 

ood e v e2 + d t v e0 + e (v eO ■ V)v e0 = -^Vn" el + e n^Vn* o 
- (E 2 + 6> e (v el x B + v e0 x Bi)), 
ud e nj| 2 + 5 t n* + 6» e V • (v el + njj v e0 ) = 0, 



ud e v i2 + d t v i0 = -a 2 Vnf + — (Ei + v i0 x B ), 



udnnL + d t nL + V • v i0 = 0. 



Because (j3.7|) implies that i^o is a gradient, the first equation implies that 

lL v e0)p 



dtBofi = 0. This yields Bq = 0, and, with (|3.6|) . we find that (rig v e o)p = 0. 
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Because v e o iP is assumed to be non identically zero (see the form of the 
initial condition (|3.3|) ). this implies n\ Q0 = 0, and finally n^ = 0. That is, 
the electronic fluctuation of density is O(e), in accordance with the rescaling 
of section 12331 For p G {-1, 1}, 

( E 2 \ = / -d t E + VxB 1 + (n« lVe0 + n» v el ) - £ Vil \ 
01 V v e2 J p \ -d t v e0 - ^(VnJi - (v e0 • V)v e0 ) - (v e0 x Bi)) J p 

In the above right-hand side, the nonlinear terms are 

(n^v e0 ) p = nl 10 v e0 , p , (v e0 x Bi) p = v e0 , p x J3 1)0 , p = -1, 1, 

and 

(v e0 • V)v e0 ) p = 0, p= -1, 1, 

a transparency relation for the convective term. The compatibilty relation 
is the Schrodinger equation for the electric field, 

1 9 

-2ipujd t E 0tP + A e £ , P - m E o,p ~ n li o e Q,p + ^ E 0p x B 10 = 0, (3.8) 

where 

A e z := 0fV(V • z) - V x (V x z). (3.9) 
The nonoscillating terms satisfy 

# e ((v e0 • V)v e0 ) = -e e Vnl 10 - (£ 2 ,o + e (v eO x Bi)„). (3.10) 

The above equality ()3.1(Jj) is the crucial equation that couples the Schrodinger 
equation (|3.8|) and the evolution equation for n\ x , which is made explicit 
below. The equations for the terms in 0(1) also contain the relation: 

1 u 

V X Bifi = Ve2,0 ~ 7-Uil,0 + (n^v e o)o, (3.11) 
Ve 

and a linear wave equation for Ujo,o an d . Because the initial data for 
Viofi and n il0 are null, n- 00 = and Vio,o = 0. With (|3.7[) . and because 
^0,0 = 0> this implies v elj o = 0. 

Equations for the terms in 0(e). 
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ud e B 3 + d t Bi + V x E 2 = 0, 

cod e E 3 + d t Ei - V x B 2 = v e3 + i%v e i + n^ v e2 + n^ 2 v e0 

- j-(yn + nljVio + nf v a ), 

udgV e3 + ^Vei + e ((v e l • V)v e0 + (v e0 ■ V)v e l) = -0 e Vn| 2 

- (E 3 + (9 e (v e2 x B + v e i x Bi + v e0 x B 2 )), 
^9 e v i3 + dtvu + (v i0 • V)v i0 = -a 2 VnJj + a 2 n| Vn| 

+ — (E 2 + x B + v i0 x Bi), 

wdflnjlg + a t n^ + e V • (ve2 + n^Veo + n^ v e i) = 0, 
w9 e nf 3 + e^n^ + V • (v a + nf v i0 ) = 0. 

In particular, 

dtvnfl + a 2 Vnf 1 = ^-£2,0, 
d tn\ lfi + V • ua,o = 0, 

<% n el,0 + ^eV • {Ve2,Q + (n* 1 V e0 )o) = 0. 

The last two equations in the above system, together with (|3,11|) . imply the 
quasineutrality relation, 

™!i,o = 4,o- (3-12) 
The first two equations in the above system, together with (|3.1U|) and H3.12|) . 
give 

J dtvn,0 + (a 2 + l)Vnf 10 = -(v e0 • V)v e0 + v e0 X Bi) , 

I d t n li,o + v ' v n,o = 0, 

where the nonlinear term can be computed with the above polarization 
conditions, 

((v e0 • V)v e0 + v e0 x Bi) Q = V|£ ,p| 2 - 
The equations at order 0(1) also yield 

d t B lj0 + V x E 2 , = 0, 

which, together with ()3.1(Jj) . implies that E20 is a gradient. Hence the term 
B\fl in (|3.8j) vanishes. 
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Finally, Eq iP and nufi satisfy the vector Zakharov system, 



1 u 

-2ipdfEo p + A e Eo p — -t^Eq p = n a q Eq_ p 

& 2 # »i,u # (313) 

(d 2 - {a 2 + l)A)n« 10 = -A|£ , P | 2 . 

In (|3,13|) . p = 1 or p = — 1. The Laplace- type operator A e was introduced 
in (j3.9j) . With the initial condition Eq : i = E, rinfi = 0, and dtnufl = 0, 
Ozawa and Tsutsumi's result ^H] guarantees the existence and uniqueness 
of a solution (Eqi, nn$) to (|3.13|) . over a time interval [0, t*), with the same 
Sobolev regularity as the initial condition. Note that the crucial coupling 
term A|£o,p| 2 comes from the convective term and from the Lorentz force 
term. In the above Schrodinger equation, the term £j-Eb,p means to a small 
shift in frequency. This term was not present in the (Z) system given in the 
introduction; it accounts for the contribution of the ions to the fundamental 
frequency. Indeed, by letting e = (cOp e io) _1 ) we took as a reference the elec- 
tronic plasma frequency u pe , which is only an approximation of the plasma 
frequency, 

47re 2 no( 1 )• 

m e rrii 

Because the ration m e /mi is equal to e 2 /0 2 (a consequence of (|1.7)) ). at first 
order in e 2 , 

(J p = Upeil + -p), 

and thus the shift in frequency in (|3.13|) means that the electric field actually 
oscillates at the plasma frequency u p . 

Higher-order terms. The WKB expansion can be carried out up to any 
order. For m > 2, the terms £ m ,p ) ?1 «,m+i,0) P = ~ lj 1 are seen to satisfy a 
linearized Zakharov system of the form 

(-2ipd t + A e - -L)E m , p = nl m+10 E , p + n\ lfl E m # + r mjP , 

"e (3-14) 

{d 2 - (a 2 + l)A)nf m+10 = -A(E miP E _ p + E 0yP E m _ p ) + r m>p , 

where r m ^ p represent the p^ harmonics of smooth functions of the profiles 
{c^Ufc/, k,k' < m— 1}. The system (|3.14j) is the linearization of (|3.13j) around 
P' n !i,o)- The hiitial data for E m ^ p and ra| rri+10 are null. Ozawa and 
Tsutsumi's method for the Zakharov equations JH] allows to solve the initial 
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value problem for ()3.14() . It provides an existence time which is a priori 
smaller than the existence time of the data r mjP , Eq iP and vtL . The solution 
has the same Sobolev regularity as the data. 

Thus, by induction, we can construct a family of profiles u m that de- 
termines an approximate solution u e app , as follows. If we assume the u^, for 
k < m — 1, to be known, with enough Sobolev regularity, then 

• E m ^ p and ni m+ \fl are defined as the unique solution of (Z) m , 

• the terms u mjP , for p G 7Z*, are deduced from E mjP and ni m+ ifi by 
polarization conditions similar to the ones found in the first terms of 
the expansion; 

• the terms u miP , for p ^ 7Z*, are computed by elliptic inversions; that 
is, in terms of (ipuj + L )~ 1 r rn _i jP . 

We obtain a profile u m , whose Sobolev regularity is smaller than the regu- 
larity of u m _i by one, and the induction is complete. The profile 

u £ = (B £ E e v e n e,fl v £ n £ 'h 

"app \ ±J ) JJ > v e ' xl e > v j ' xl i 7 ) 

provides, via 1)3.4(1 . an approximate solution to (EM)". Then 

< = [B £ , E e , v £ , i log(l + en^), vf , ~ log(l + en^% =w(/e2 , 

is an approximate solution to ()3.2() . in the sense of ()2.19() . The approxi- 
mate solution can be made arbitrarily precise, in a Sobolev norm, provided 
the initial data have enough Sobolev regularity. Then, by construction, u £ a 
satisfies all the conditions stated in Assumption 12 . 1 fl 

3.2 Stability of the approximate solution 

Consider now a perturbation of the initial condition a introduced in ()3.3() : 

o £ :=a + £ fc V, (3.15) 

where ip e is a bounded family in H°~ l °~ 2 , and ko + 3 < Iq < a — 2 — |. We 
need to assume that (jl.4j) is satisfied so that a £ is a proper initial datum for 
the Euler-Maxwell system. This amounts to the conditions, 

V • B^e =0, V • E^e = — (n eilp e - U^), 

for the coordinates of the perturbation ip £ . 

Consider the approximate solution u £ a to 1)3.2(1 , associated with the initial 
datum a, at order Iq. 
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Theorem 3.1. If ko > 3 + |, the system (|3.2|) . together with the initial 
datum (|3.15j) . has a unique solution u £ defined over [0, t*), independent of 
e. Moreover, for all < to < t* , there holds, for e small enough, 

sup \\u £ - ul\\ £>s < Ce^ 0-1 , 
o<t<t 

with a constant C independent of e, and a Sobolev index s > f • 

The proof shows that the existence time t* is bounded from below by 
the existence time of the approximate solution u £ a . Note that the estimate 
()3.16j) . the condition > 3 + |, and the description of the approximate 
solution given in the above section, imply 

sup sup (\E £ - (E 0il e iujt/£ ' 2 + c.c.)\ + \n £ - en ilfi \) < Ce 2 , (3.17) 

0<t<io xeR 3 

where E £ , n £ are coordinates of the solution u £ of (|3.2j) , as in 1)3. lj) , where 
n H,0 = \ log(l + enJ 1 ), and I?o,i and -n\ x solve the Zakharov system (|3,13(1 . 

with the initial condition Eo,i(t = 0) = E°, n\ lQ (t = 0) = 0, dtn\ lQ (t = 
0) = 0. The asymptotic estimate (|3.17|) is thus the estimate that validates 
the Zakharov model, as it actually gives a description of the electric field 
E £ and the fluctuation of density n £ in (EM) by means of the solution 

(£0,1,^1,0) of ( Z V 

Theorem 13 . 1 1 follows as a corollary of Theorem 12. 151 if one can prove that 
the (EM) system satisfies the assumptions of Theorem l2.15l An approximate 
solution satisfying Assumption 12.111 was constructed in the above section. 
The next sections are devoted to the verification of Assumptions I2.1UI (hy- 
perbolic structure, regularity of the eigenvalues and eigenproiectors). 12.121 
(localization of the resonances) and !2.13l (transparency). 

3.2.1 Eigenvalues and eigenvectors 

We check in this section that the (EM) system satisfies Assumption I2.1UI 
The operator A defined at the beginning of section El obviously belongs 
to C^M 1 . For all e,u,^, A(e,u,^) is a hermitian matrix, and satisfies a 
decomposition of the form (|2.16|) . It remains to check hypothesis (i) to (iii) 
in Assumption 12.1(71 

The eigenvalues arising in the spectral decomposition (|2.16|) are the so- 
lutions, 

qj = A(s,u,€), e > 0,u € R 14 ,£ G M 3 , 



(3.16) 
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of the polynomial equation in u, 

det (iw + A(e, u, = 0. (3.18) 

A look at the definition of A\ shows that the eigenvalues depend on u only 
through the scalar terms, 

x := e6 e v e • £, and y := e 2 Vi • £, 

representing the electronic and ionic convections (v e and Vi are coordinates 
of u). Equation (|3.18j) factorizes into a transverse, degree four equation, 

2 2 

(u> - x)(u, - y)(^ 2 - 1 - |£[ 2 - £ ¥ ) = x(u; - y) + ^y(w - x), (3.19) 
and a longitudinal, degree five equation, 

W (( W - y) 2 - . 2 a 2 |e| 2 )((^ - x) 2 - 1 - fl e 2 |£| 2 ) = -x(( w - y) 2 - £ 2 a 2 |£| 2 ) 

+ ^(u,-y)(( w -x) 2 -0 2 |£| 2 ). 

(3.20) 

For all e, x, y, £, the Kernel of w4(e, it, £) has dimension one. It is generated 

by 

e := (Jy, 0,0, 0,0,0). 

The solutions of (|3.19j) have (algebraic and geometric) multiplicity two, while 
the solutions of (|3.20f) have multiplicity (algebraic and geometric) one. The 
solutions of 1)3. 19j) and (|3.20j) are algebraic functions of e, £, x, y. Evaluations 
of these functions uj = u>(e, (x, y),£) at (x, y) = (e6 e v e ■ £,£ 2 Vi ■ £) give the 
non-zero eigenvalues of the (EM) system. 

For (x,y) = (0,0), equation (|3~T9~]) and (|3~2l7|) simplify to 

cu 2 (a, 2 -l-|e| 2 -i=0, (3.21) 

and 

- e 2 a 2 |e| 2 )(- 2 - 1 - 9l\tf) - - 9l\tf)) = 0. (3.22) 

The solutions u> = u>(e,£) of (|3.21[) - l)3.22|) are represented on figured 
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Figure 2: The characteristic variety for the operator linearized around 0. 
The transverse modes at (e, 0,£) are 



A s (e_) = 0, A s (e+) = 0, A± = ±^ 1 + |£|2 + -. (3.23) 
The Klein-Gordon, longitudinal modes at (e, 0,£) are 

, ± = ± v / rrw + 2( , g(1+ e ^ |2) + o^), 

and the acoustic, longitudinal modes at (e, 0,£) are 

fi a (e-) = 0, M s (e+)± = 0(e), ^ = ±^1 + 6^ + 0(6), 

locally uniformly in £. The Klein-Gordon modes ji± have constant multiplic- 
ity, hence are analytical in e, £. Crossing of eigenvalues occur for the acoustic 
modes only at the zero frequency Thus away from the zero frequency, the 
acoustic modes are analytical as well. It is easy to check that 

0<Me + )±|<Ce|£|, (3.24) 
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uniformly in e, £, and that 



(e+)± = ±e|e|(Vl + « 2 + 0(e 2 + |£| 2 )), (3.25) 



uniformly in e, |£| < cj. Regularity at infinity can be directly checked using 
an exact description of the solutions of (|3.22|) : one finds that the longitudinal 
eigenvalues have the form 

1 1 £ 1 1 £ 

M s (e+)± = ±ea|f| + j^|-F±(|^|,-j^j-,e), M± = e |£| + j£j-G±(-j|p j^pe), 

(3.26) 

where F± and G± are analytical in (|£o|,oo) x S 2 x R + , for |£q| > 0. 

Thus the eigenvalues at (x, y) = (0, 0) satisfy (|2.1j) (regularity in e, £ for 
£ away from zero, decay at infinity) with m = 1, that is, condition (ii) in 
Assumption 12 . lUl is satisfied by the eigenvalues. 

For £ in a compact subset of M 3 , the eigenvalues evaluated at (x,y) = 
(s6 e v e -t;,E 2 Vi-£,) are small perturbations of the eigenvalues at (x, y) = (0,0). 
Thus the Klein-Gordon modes are separated from the acoustic modes. For 
large £, the contribution of the convective terms to the eigenvalues is not 
negligible, but the acoustic modes are all 0(e|£|), while the Klein-Gordon 
modes are 0([£|). Thus condition (i) in Assumption 12 . 1 ( )l is satisfied. 

The Klein-Gordon modes are single eigenvalues of ()3.19|) and (|3.20[) . for 
all £. Hence the eigenvalues and the eigenprojectors corresponding to the 
Klein-Gordon modes are analytical 

For £ 7^ 0, we let {£1,^2} be an orthonormal basis of {^} ± . 

At (e, 0,£), £ 7^ 0, the eigenvectors associated with the transverse eigen- 
values are 

6± - ( ^' ei '^' '^A^' 0) ' 



e 



1,1 A± ' s ^ ' A± ' e A± 

e s (e") := -==L= (£i,0,»|£|f 2 , 0,0,0) : 

v 1 + l£r 

e' s (e-) := J_ (g 2 , 0, 0, 0, 0) ; 

V 1 + Iff 

e - (e+) := rari ( W' ' ' ' ft ' 0) ' 
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At (e, 0, £), £ 7^ 0, the eigenvectors associated with the longitudinal eigen- 
values are 



/±:= 



and 



/4~ W £ -^e|£| ie jj^ff £ M± 

iae 2 jgj 
# e ii± — e 2 a 2 \r]\ 2 fi± 

f s (e~) := (0, -tfl e g, 0, 1, 0, ), 



, , +, •< - , +. £ -i/9 e /2 s (e+) ± /z s (e+) ± £ 



e M-C^-flgiei 2 )^ ' 1*1 ' ^(e + )±^' 
where 

- , + , ^(e+^-e 2 ^ 2 

From (|3.25|) . one sees that ji s {e + )± = 0(e|£|). This yields a description of 
/s( e+ )± f° r small frequencies: 

/ s (e + ) ± = (0, 0(|£|), O(e), 0(1), O(l), O(l)) . 

The corresponding orthogonal eigenprojectors (TIj , 1 < j < n with the no- 
tation introduced in Assumption 12,10^1 are defined by 

e ® e + e' <g> e', / <g> /, 

where e = e±,e s (e~) or e s (e+), / = f±,f s {e~) or / s (e+)±, ande,/ := jfy, rjj- 
With the above description of the eigenvalues for small and large |£|, it is 
straightforward to check that all the eigenprojectors, evaluated at (e, 0), are 
Fourier multipliers in the class C°°.M . Thus condition (ii) in Assumption 
12. 101 is satisfied. 

Condition (iii) follows from the above description of the acoustic modes, 
as they all have a prefactor e. 

We now turn to condition (|2.18j) . The total eigenprojectors, defined in 
(t2~T71) . are 

n = S± ® £t + S± ®e± + / ± ® /±> 
± 
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and 

n s = X^( e± ) ® e s (e ± ) + e^(e ± ) ® e' s (e d 
± 



+ ® />") + e ® e + ^ / s (e+)± ® / s (e+)±. 

± 

Given (eo,u ,{;o) G [0,1] x M 14 x M 3 , for (e, u, £) in a neighborhood of 
(eoi^QjCo), 

n (e,u,£) = — f (z-A(e,u,())- 1 dz + ^- I (z - A(s, u, O) -1 ^, 

(3.27) 

and 

n.(e, «, = ^ / - «, (3-28) 

where To+ (resp. To-) is a contour enclosing all the positive (resp. negative) 
Klein-Gordon eigenvalues at (eo,uo,^o), and no other eigenvalues of A, and 
T s is a contour enclosing all the acoustic eigenvalues at (eo,^o>£o)> an d no 
other eigenvalues of A. The above description of the eigenvalues show in 
particular that such contours do exist. It follows from these representations 
that Ho and Il s are analytical in e, u, £. To investigate the behaviour for 
large £, notice that A can be written, for |£| ^ 0, 

A(e,u,£) = Aq{e) +A 1 (e,u,0, 
where A\ is linear in £, so that, using polar coordinates £ = 9) £ x§ 2 , 

A(e,u,0 = \Z\(jzAo(e) + A 1 (s,u,9)). 

Thus an eigenvector e(e, u, |£|, 0) of A(e, u, £), associated with an eigenvalue 
X(e,u, equals an eigenvector e(e,u, t|t>#) of 

A(e, u, |£|, 0) := lei^Ce) + A l (e, u, 9), (3.29) 

associated with the eigenvalue ^X(e,u, ■mjQ)- The operator A is the long- 
wave operator associated with A, introduced in [27j in the study of the 
short-wave limit. 

For small |£|, the eigenvalues of A split into "acoustic" eigenvalues, of 
size 0(s), and "Klein-Gordon" eigenvalues, of size O(l). The total eigenpro- 
jectors, 

n = -L / (z - A(e, u, \^\,9))- 1 dz + J- / (z - ^(e, u, fl))" 1 ^, 
2*w 7f o+ 2 i vr 7 ro _ 
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and 

are analytical in e, u, for |£| in a neighborhood of 0. There holds 

1 £ 

IIj(£,u,£) = n 3 -(e,w, , ), j = 0,8, 

and thus the total eigenprojectors IIo,n s of A are analytical in e,u, 4r, t|j, 
for £ in a neighborhood of oo. This implies that they satisfy decay estimates 
of the form (|2.1|) . and (j2.18j) is proved. 

To conclude this section, we now indicate how the Klein-Gordon eigen- 
values and eigenvectors depend on x and y. These descriptions are needed 
in the evaluation of the interaction coefficients that enter Assumption 12.131 

The eigenvalues satisfy 

d * x ± = ^rrw) +0{£% ^ ± = 27irtm + ° (e) ' 

at (x,y) = (0,0), locally uniformly in £. 
The eigenvectors are 

f\v\V2 -irn ie m n \ 

4 == (^^,^,0,^^,0), 



and 



X± X± — x 9 e X± 



is (//± - x) 2 - 9 2 \i]\ 2 r] /j,±-y 
9e (/u± - y) 2 - £ 2 a 2 \r]\ 2 \rj\ fj,±-x' 
ie 2 (n± - x) 2 - 6 2 \r}\ 2 \rj\ 



(fx± — y) 2 — e 2 a 2 \rj\ 2 [j,± — x7 



3.2.2 Resonances 

We use here the description of the eigenvalues given in section 13.2.11 The 
resonance equation Q2.2UJI for j,k < no (Klein-Gordon/Klein-Gordon reso- 
nances) and e = is 

vTTf-ViTW-uo. (3.30) 
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— h- 

0(1) 



Figure 3: (0-0) resonances. 



If 9 e is small enough, we choose c m = 1 and C m = 2. The left-hand side 
in (|3.30|) is bounded away from zero for |£| ^ [1,2], and hypothesis (0-0) in 
Assumption 12. 121 is satisfied. 

The resonance equation (|2.20|) for j < uq < k (Klein-Gordon/acoustic 
resonance) is 

V/1 + K 2 |e 2 |-1-M = 0. (3.31) 

where k = 1 or k = 9 e , and fj, = or fi = ^ s {^ + )±- If &e and £ are small 
enough, we can choose q = 1/2. Then the left-hand side in ()3.31|) is bounded 
away from zero for |£| > 1/2. 

The resonance equation (|2.21|) (secondary Klein-Gordon/acoustic reso- 
nances) is 

\A + « 2 I£I 2 -(P + P') = 0, (3.32) 

where p,p S {—1,1}, k = 1 or k = 8 e . The left-hand side in (|3.32|) is 
bounded away from for |£| < 1, and hypothesis (0-0-s) in Assumption ^. 121 
is satisfied. 

The resonances are pictured on Figures to |SJ 
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Figure 4: (0-s) resonances. 



3.2.3 Transparency 

We check in this section that Assumptions 12. l"3l and U. 141 are satisfied. 

With the definition of A given at the beginning of section |HJ one has 
A^ = 0. It is straightforward to check that II s (0, 0)V(u a )n (0, 0), where V 
is defined by Q2.23J1 . has size 0(e), and Q2.27JI is satisfied. 

Up to terms of size 0(e 2 ), the source term defined in (|2.22j) is 



It depends on the coordinates B, v e , n e , of the approximate solution u a , 
evaluated at 9 = cot/e 2 . 
One computes, 





n e v e + en el v e 
6» e (v e x B + (v e • ()v e + eBi x v e + ev e ■ Vv e ) 
-^ e (v e • £,)n e 



\ 



B(u a )u 



\ 








/ 



e* + B(u a )e s (e ) 




+ 0(e) 



64 




and also 

f+B(u a )fs(e ) = 2 \7T + 

Is I 

In particular, these interaction coefficients are not small for small frequen- 
cies. This shows that the non-transparency condition (|2.24j) is satisfied. 

We now check that the transparency condition of Assumption 12.131 is 
satisfied. 

The symbol p defined in Q2.25JI is p = p^ + 0(e), where 

P w : =n s (^(n s i (0) )9 t n -^ a + ^n s /)(9> a ) 

|a|=l 

+ u s A^d^u d v u ■ d a x u a - d?n s B{o, 9> a ))n , 

where, unless otherwise noted, the symbols and their derivatives are evalu- 
ated at (e, v) = (0, 0). 

Assumption 12.131 is a transparency assumption for the interaction co- 
efficient d u B r (e, 0) • u a . Direct computations, using the description of the 
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eigenvalues and the eigenvectors given above, yield the bound 

|n fc (e,0)(a nJ B r (e,0)- Ua )n j ( £ ,0)| <C B (|e| 2 + e|el), (3.33) 

for j < tiq < k, uniformly in e,t,x and |£| < q. Now (|3.33[) and the above 
description of the phases imply (|2.28|) . as follows. Let 5 > be given. 

• If |f | < 5e, then (EHHl) directly implies (l2~m with C = C B S(1 + 5). 

• If |f| > Se, then for < e < e , |*j,fc )P (e)| > Co(<5)|£| 2 , with 

3fl e 2 C(a,e ) 
c o(<3) := — 



6 ' 

for some nondecreasing function C(a,eo). If ^ is chosen to be large 
enough, then Co (5) > 0, and (|2.28|) is satisfied with 

C = C B (- X + 1)C (5)- 1 . 
o 

Assumption 12.141 is a symmetrizability condition for the interaction co- 
efficients 

B = U B(u a )U , B s = U S BU S . 

As these interaction coefficients enter the equation with a prefactor ^, it 
is sufficient to consider the leading term in e in B(u a ). In particular, the 
contribution of in Bq and B s is only 0(e). We can write indeed 

op 6 (^ 1 )(u))w« = eo Pl (A {1) (u))u a , (3.34) 

because is linear in f , and the norm of 1)3. 34 j) is bounded by 

(\\u a \\i )S+do )\\u\\ £ , s . 

The auto-interaction coefficients are all purely imaginary: 

e*B{u a )e e iR, 

where e is any eigenvector of the Euler-Maxwell equations. The other inter- 
action coefficients between the Klein-Gordon modes are 

e*±B(u a )e'± = 0, e* ± B(u a )f± = -id e ^ e ■ 

I' : 

(e' ± )*B(u a )e ± = 0, (e' ± )*B{u a )f± = -i9 e ^v e ■ 

f* ± B{u a )e ± = -i6 e ^v e ■ 6, r ± B{u a )e' + = -i6 e ^v e • 
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where for all £ 7^ 0, {^1,^2} is an orthonormal basis of {£} . Let So = 
5 o +) + < s o~ ) > where 

S { ±] := e± <g> e± + e' ± <g> e± + — /± ® /±. 
Then So is a symmetrizer, in the sense that 

i(S 5 + (SoA))*) G C°°M°. 

Finally, we turn to the interaction coefficients between acoustic modes. Re- 
mark first that the divergence-free condition for the magnetic field, 

V • B = 0, (3.35) 

is equivalent to 

op e (e )u = 0. 

Condition (|3.35|) is propagated by the equations, that is, the solution belongs 
to the orthogonal of the image of op £ (eo) in L 2 (R 3 ) if the initial datum does. 
Thus we can overlook the interaction coefficient involving eo- At first order 
in e, the interaction coefficients with the other acoustic modes all vanish, 
except for 

Ue„yB( Ua )f s (e +)± ^ e+)± -^" 2 



(/ S (e + ) ± )*£K)/ S (e_) 



VTTW ^ 



It follows from the description of the eigenvalues given in sect ion ?S . 2 1 1 1 that . 
for e small enough, 

1 < -0m 2 < 2 

2- / i,(e+)|-^|2- • 
uniformly in £ £ IR 3 . Thus we can take 7 = 2 in Assumption 12 . 141 Let 

S s :=e <8> e + e s (e_) <g) e s (e_) + e s (e + ) ® e s (e+) + / s (e_) <g> / s (e_) 

02ICI2 



□2 1 1 |2 



+ ■- ~ ^ £ _ fl2 |.,| 2 ^( e +)+ ® /«( e +) + 
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Then 



-(S s B s + (S s B s )*)eC co M°. 



Because Sq and S s are diagonal in a basis of eigenvectors of A, the matrices 
SqAq and S S A S are hermitian. Finally, 



defines a symmetrizer that satisfies Assumption 12. 13l 
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